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5_^ ■ We present gas mass fractions of 38 massive galaxy clusters spanning red- 

O^- shifts from 0.14 to 0.89, derived from Chandra X-ray data and OVRO/BIMA 

^ ' interferometric Sunyaev-Zel'dovich Effect (SZE) measurements. We use three 

models for the gas distribution: (1) an isothermal /3- model fit jointly to the 

X-ray data at radii beyond 100 kpc and to all of the SZE data, (2) a non- 

Q>^ ' isothermal double /3-model fit jointly to all of the X-ray and SZE data, and 

^ ■ (3) an isothermal /?-model fit only to the SZE spatial data. We show that the 

Tj- ■ simple isothermal model well characterizes the intracluster medium (ICM) out- 



9"' 
6 



cd 



side of the cluster core, and provides consistently good fits to clusters spanning 
O ■ a wide range of morphological properties. The X-ray and SZE determinations 

of mean gas mass fractions for the 100 kpc-cut isothermal /?- model are /gas(X- 
ray)= 0.110tg:[iE]i t[]:E]?^ and /gas(SZE)= OA16t°oml-oml, where uncertainties are 
statistical followed by systematic at 68% confidence. For the non-isothermal dou- 
■^ . ble /3-model, /gas(X-ray)= 0.119+°;°°^ +°:°?I and /gas(SZE)= 0.121+^05 -oml- For 

the SZE-only model, /gas(SZE)= 0.120t[|:|][|^ tg;!]^?. The agreement in the results 
shows that the core can be satisfactorily accounted for by either excluding the 
core in fits to the X-ray data (the 100 kpc-cut model) or modeling the intracluster 
gas with a non-isothermal double-/? model. We find that the SZE is largely insen- 
sitive to structure in the core. Our results indicate that the ratio of the gas mass 
fraction within r25oo to the cosmic baryon fraction, q T\ — , is 0.68;!^Q;]^g where 
the range includes statistical and systematic uncertainties at 68% confidence. 
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Finally, by assuming that cluster gas mass fractions are independent of redshift, 
we find that the results are in agreement with standard ACDM cosmology and 
are inconsistent with a flat matter dominated {Qm = 1) universe. 

Subject headings: cosmology: observations; cosmological parameters; X-rays: 
galaxies: clusters 



1. Introduction 

Massive galaxy clusters (M ~ 10^^ M©) are thought to be the relatively recent descen- 
dants of rare high density fluctuations in the primordial universe. The evolution of massive 
clusters is critically dependent on cosmology, in particular on the matter density Qm and 
the normalization of the power spectrum erg. Furthermore, as massive galaxy clusters have 
collapsed from volumes of order 1000 Mpc^, their composition should reflect that of the non- 
relativistic components of the universe, i.e., the baryon budget of a cluster should reflect 

For these reasons there is considerable interest in using cluster observations to constrain 
cosmology, and there has been significant success. For example, early X-ray observations 
were used to constrain the gas mass fraction of clusters and thereby set a lower limit to 
Qb/^m- Using constraints on Qb from Big Bang Nucleosynthesis, this led to upper limits 
on Qm that strongly ruled out a flat matter dominated universe (White et al. 1993); a precise 
measurement of Qm was not possible due to insufficient understanding of the cluster baryon 
budget. A low Qm was also indicated by the discovery of high redshift (z ~ 1) massive 
clusters, which should be extremely rare in a matter dominated, flat universe (Bahcall & 
Fan 1998; Donahue et al. 1998). 

Cosmological constraints can be obtained by exploiting the expected redshift indepen- 
dence of cluster gas mass fractions; only if the correct cosmology is used in the derivation of 
the gas mass fractions for a sample of clusters spanning a reasonable redshift range would 
the resulting gas mass fractions be constant with redshift (Sasaki 1996; Pen 1997). This 
technique is independent of the uncertainty in the cluster baryon budget as well as the value 
of the Hubble constant. 

Subsequently there have been many more cosmological studies using cluster X-ray and 
Sunyaev-Zel'dovich Effect (SZE) measurements (e.g., see White et al. 1993; David et al. 
1995; Evrard 1997; Myers et al. 1997; Mohr et al. 1999; Ettori & Fabian 1999; Grego et al. 
2001; Allen et al. 2002; Sanderson & Ponman 2003; Allen et al. 2004). In general the results 
are in good agreement, along with cosmic microwave background (CMB) (i.e., Spergel et 
al. 2006) and Type la supernova measurements (Perlmutter et al. 1999; Riess et al. 1998), 
with the now standard ACDM cosmology. 



Future large scale cluster surveys are being pursued which exploit the critical cosmo- 
logical dependence of the evolution of cluster abundance to provide independent, precise 
determinations of Qm, ^a, o"8 and the equation of state of dark energy, w{z) (e.g., Kneissl 
et al. 2001; Romer et al. 2001; Kosowsky 2003; Ruhl et al. 2004). The ability to extract 
cosmology from large surveys using SZE, X-ray, or any other cluster observable, will depend 
on the ability to link the measurements to cluster mass. An improved understanding of 
cluster structure and evolution is of interest in itself and would clearly also be helpful for 
improving cosmological measurements with clusters. 

In this paper we investigate the gas mass fractions for a sample of 38 clusters spanning 
redshifts from 0.14 to 0.89, using Chandra X-ray data and BIMA/OVRO interferometric 
radio SZE data. This SZE cluster sample is the largest yet compiled by over a factor of two. 
Combining the SZE data with the high resolution Chandra X-ray data we are able to compare 
the SZE and X-ray gas mass fractions using two different models for the ICM. Comparison 
of the SZE and X-ray results allows us to place constraints on the possible systematic effects 
associated with our incomplete knowledge of cluster properties. The results also provide 
constraints on the cluster baryon budget, and thereby provide an important observational 
benchmark for cluster simulations. The baryon budget is sensitive to the details of cluster 
formation, gas cooling and star formation history of the ICM (Ettori et al. 2006). Lastly, we 
consider the cosmological constraints derived by assuming the cluster gas mass fractions of 
our sample do not evolve with redshift. 

The paper is organized as follows: we review the theory underlying the X-ray emission 
and the SZE in clusters in §2 and describe the X-ray and SZE data in §3. In §4 we describe 
the cluster models and analysis methods used to determine the gas mass fractions, and we 
present tests of these models. In §5 we discuss additional sources of uncertainty in the gas 
mass fraction measurements, both statistical and systematic. The results are presented in §6, 
and constraints on the cluster baryon budget and cosmological parameters are detailed in §7. 
Finally, we summarize our conclusions in §8. All uncertainties are at the 68.3% confidence 
level, and a ACDM cosmology with ^2^=0. 3, Qa=0.7, h=0.7 is assumed unless otherwise 
stated. 



2. X-ray emission and Sunyaev Zel'dovich Effect 

X-ray emission in clusters arises predominantly from thermal bremsstrahlung for gas 
electron temperatures Te > 3 keV. The X-ray emissivity is not a directly observable quantity; 
instead X-ray proportional counters and CCDs measure the X-ray surface brightness over 
some frequency band (e.g., Sarazin 1988), 

^^ = A /I , — \d / neUuAeH d£ (1) 

47r(l + z)*J 



where the integral is along the line of sight and A is called the X-ray cooling function and 
is proportional to Tg . The observed X-ray emission can also be used to measure the gas 
temperature using instruments with spectral capability. 

The thermal SZE is a small (< 1 mK) distortion in the CMB spectrum caused by inverse 
Compton scattering of the CMB photons off of energetic electrons in the hot intracluster gas 
(Sunyaev & Zel'dovich 1970, 1972). The spectral distortion can be expressed for dimension- 
less frequency x = hv/kjiTc^i^ as a temperature change AT relative to the CMB temperature 

J- CMB- 

AT „_ „_ /• KT, 



rj. f{x)y = f{x) a^Ue ^de, (2) 

where y is the Compton-y parameter, a^ is the Thomson scattering cross-section of the 
electron, f{x) contains the frequency dependence of the SZE, and the integral is along the 
line of sight. The frequency dependence can be expressed as 

/(^) = (^ J^T - 4) (1 + 5szE(a:, Te)) (3) 

where 6szE{x,Te) is a relativistic correction. We adopt the analytic corrections of Itoh et al. 
(1998) which are good to fifth order in k-sTe/meC^. The linear density and temperature 
dependences of the SZE make it a complementary probe to the X-ray emission, which varies 
as UgTe . The SZE is a decrement at low frequencies (< 218 GHz) and an increment at 
high frequencies due to the up-scattering of photons by the hot electrons. There is also a 
kinetic SZE that arises from scattering of CMB photons in a cluster with line of sight motion 
relative to the CMB rest frame. Discussion of the SZE in this work refers to the thermal 
SZE unless otherwise specified. We address the possible effects of the kinetic SZE on our 
results in §5. 

With knowledge of the gas temperature, one can use either Sx or AT to determine the 
gas distribution and obtain Mgas. The total mass Mtotai can be determined by assuming that 
the gas is in hydrostatic equilibrium (HSE) with the cluster potential. The gas mass fraction 

/gas is then Mg^s /Mtotai- 



3. Data 

The cluster sample consists of clusters that have both X-ray data from the Chandra 
X-ray Observatory and SZE data from the BIMA/OVRO SZE imaging project, which uses 
the Berkeley-Illinois-Maryland Association (BIMA) and Owens Valley Radio Observatory 
(OVRO) interferometers to image the SZE. The BIMA/OVRO SZE imaging project has 
generally pursued hot clusters (published Tx ^ 5 keV) which by inference should be mas- 
sive and have a strong SZE. Properties of the cluster sample, including full cluster names, 
redshifts, positions, and Chandra and BIMA/OVRO observing information, are listed in 
Table 1. 



The X-ray data were reduced using the Chandra Interactive Analysis of Observations 
(CIAO) 3.2 software package with CALDB 3.1. Details of the data reduction procedure 
are given in Bonamente et al. (2004) and Bonamente et al. (2006); briefly, CIAO and the 
CALDB are used to apply corrections for charge transfer inefficiency (CTI) and buildup of 
a contaminant on the optical blocking filter. The observations are filtered to contain only 
photons with energies between 0.7 and 7.0 keV, and to remove contamination from solar 
fiares. The data are then binned into images for each chip relevant to the observation, and X- 
ray backgrounds are extracted from appropriate cluster-free regions of the chips (Bonamente 
et al. 2004). Exposure maps are also created with CIAO, and all spatial analysis is performed 
on the binned images and exposure maps. Cluster spectra are extracted from a circular region 
centered on the cluster containing roughly 95% of the cluster counts, again restricted to the 
energy range 0.7 — 7.0 keV. Spectra are grouped so that Gaussian statistics can be used in 
the spectral fitting, and a background spectrum is also extracted using the same background 
regions used in the image analysis. The software package XSPEC (Arnaud 1996) is used to fit 
a model to the spectrum and determine the X-ray spectroscopic temperature and metallicity 
of the cluster. 

Interferometric radio observations of the cluster SZE were performed at the Berkeley- 
Illinois-Maryland Association observatory (BIMA) and at the Owens Valley Radio Obser- 
vatory (OVRO). The millimeter- wave arrays were equipped with 26-36 GHz receivers for 
the SZE observations (Carlstrom et al. 1996). Most of the OVRO and BIMA telescopes 
were placed in a compact configuration to provide sensitivity on angular scales subtended 
by distant clusters (typically ~ 1') and a few telescopes were placed at longer baselines for 
simultaneous point source imaging (Reese et al. 2002). 

The SZE data consist of the position in the Fourier domain [u-v plane) and the visibilities 
— the real and imaginary Fourier component pairs as functions of u and v, which are the 
Fourier conjugate variables to right ascension and declination. The effective resolution of 
the interferometer, the synthesized beam, depends on the u-v coverage and is therefore a 
function of the array configuration and source position. A typical size for the synthesized 
beam for the short baseline data is ~ 1'. Details of the SZE data reduction can be found in 
Grego et al. (2001) and Reese et al. (2002). Briefiy, the SZE data were reduced using the 
MIRIAD (Sault et al. 1995) and MMA (Scoville et al. 1993) software packages. Absolute 
flux calibration was performed using Mars observations adopting the brightness temperature 
from the Rudy (1987) Mars model. The gain was monitored with observations of bright radio 
point sources, and remained stable at the 1% level over a period of months. Data were excised 
when one telescope was shadowed by another, when cluster observations were not bracketed 
by two phase calibrators, when there were anomalous changes in the instrumental response 
between calibrator observations, or when there was spurious correlation. 

Images were made with the DIFMAP software package (Pearson et al. 1994) to inspect 
the data quality and, using only long baseline data, to identify and fix the positions of radio 



point sources. The point source fluxes are included as free parameters in the model fitting, 
using the same methodology as Reese et al. (2002). 



4. Analysis Methods and Modeling 

4.1. Cluster Density Models 

The isothermal /9-model has frequently been used in the analysis of X-ray and SZE 
cluster images (Cavaliere & Fusco-Femiano 1976, 1978; Jones & Forman 1984; Elbaz et al. 
1995; Grego et al. 2001; Reese et al. 2002; Ettori et al. 2004). The 3-dimensional electron 
number density is given by 

/ ^2 \ -3/3/2 

ne(r) = ne„ U- + ^] , (4) 

where Ue is the electron number density, r is the radius from the center of the cluster, Vc is 
the core radius of the intracluster medium (ICM), and /9 is a power law index. A convenient 
feature of the isothermal /3-model is that the X-ray surface brightness and SZE decrement 
profiles take simple analytic forms 

/ ^2\{l-6/3)/2 

S. = S.o[l + ^j , (5) 

n2\ (l-3/3)/2 



AT = ATo(^l + -J , (6) 

where Sxo is the central X-ray surface brightness, ATq is the central thermodynamic SZE 
temperature decrement /increment, and 6c is the angular core radius of the cluster. 

However, recent deep X-ray Chandra observations and numerical simulations indicate 
that the /3-model is not a good description in the outskirts (r > 1 — 1.5r25oo) of clusters 
(Borgani et al. 2004; Vikhlinin et al. 2005); to avoid biases associated with this effect, we 
compute masses enclosed within r2500 5 the radius at which the mean enclosed mass density 
is equal to 2500pcrit- Results are not extrapolated beyond this radius. Incidentally, r25oo 
is also the outer limiting radius at which both the Chandra (e.g., Allen et al. 2004) and 
BIMA/OVRO data (Grego et al. 2001) provide strong constraints on the ICM model. 

In some clusters the isothermal f3 model fails to provide a good description of the X- 
ray surface brightness observed in the cluster core. This is the case, for instance, in highly 
relaxed clusters with sharply peaked central X-ray emission. We have therefore developed 
two extensions of the isothermal /3-model to overcome this limitation; we describe these new 
models and their application to the X-ray and SZE data below. 



4.1.1. The 100 kpc- Cut Model 

First we consider a single isothermal /9-model fit to the X-ray data with the central 
100 kpc excised. The 100 kpc radius is a good compromise, as it is large enough to exclude 
the cooling region in cool-core clusters while keeping a sufficient number of photons to enable 
the mass modeling. The X-ray spectroscopic temperature is also determined using photons 
extracted from a radial shell between 100 kpc and r25oo- This is referred to as the 100 kpc-cut 
model in the remainder of this work. 

There is no simple way to excise the central 100 kpc from the interferometric SZE data, 
because these data are fit in the Fourier plane. We therefore fit the entire SZE dataset, while 
using the X-ray spectroscopic temperature from the 100 kpc-cut model. The inclusion of 
the dense core in the SZE data should have little effect on the derived cluster parameters, 
because the SZE as a probe of pressure is less sensitive to behavior in the dense core than 
are the X-ray data. It should also have little effect on the best-fit shape parameters, 6c and 
P, because these fits are driven mainly by the X-ray data. The SZE data therefore mainly 
constrain the overall normalization of the SZE signal, which is insensitive to the details of 
the cluster core (e.g., Nagai 2005). 



4.1.2. The Non-isothermal Double (3 -Model 

We also develop a more sophisticated cluster plasma model that takes into account 
temperature profiles. A motivation for considering this model is to assess the biases arising 
from the isothermal assumption and the effects of the core exclusion in the previous model. 

The model uses a second /3-model component to describe the sharply peaked X-ray 
emission present in the cores of some clusters (Mohr et al. 1999). The 3D density profile of 
the double-/? model is expressed by. 



riAY] = n^. 



^2 \ -3/3/2 . 2 \ -3/3/2- 

/(1 + - +(1-/) 1 + - 



(7) 



where the two core radii, r^ and Vc^, describe the narrow, peaked central density component 
and the broad, shallow outer density profile, respectively, and / represents the fractional 
contribution of the narrow, peaked component to the central density ne„ (0 < / < 1). This 
model has enough freedom to simultaneously fit the X-ray surface brightness in the outer 
regions and the central emission excess seen in some clusters. We set / = (equivalent of 
the single /5-model) if the reduced x^ for such a fit is less than 1.5. 

The 3D temperature profile is modeled assuming that the ICM is in hydrostatic equi- 
librium with a NEW dark matter density distribution (Navarro et al. 1996, 1997), 

PDM(r) ^ 6c , . 

Pcnt(^) (r/r,)(l + r/r,)2 ^^ 



where Pcrit(-2) is the critical density of the universe at redshift z, r^ = r2oo/c is the scale 
radius, c is the concentration parameter of the dark matter, and 5c is the characteristic over- 
density of the halo and is related to the concentration parameter c. Using the best-fit 3D gas 
density model, we solve the hydrostatic equation iteratively for the 3D model temperature 
profile. This 3D temperature profile is then weighted by the X-ray cooling function and 
the square of the cluster density, convolved with the appropriate instrumental response, and 
integrated along the line of sight, yielding a projected temperature profile model. Finally, 
the projected temperature profile model is compared to the observed temperature profile 
extracted from the Chandra data. 

Altogether, the model is parametrized by the combination of the ICM density model, 
Uco, ^ci, ^C2) l^i and /; the concentration parameter of the dark matter, c; and the outer 
radius, r2oo- This model is referred to as the non- isothermal /3-model in the remainder of 
this work. 



4.1.3. The SZE- only Model 

We also adopt a third model, used to fit the SZE spatial data independently of the 
X-ray spatial data. We use the full SZE dataset with no 100 kpc exclusion, but use the 
100 kpc-cut X-ray spectral data to determine the gas mass fractions. The SZE spatial data 
alone do not have sufficient resolution to provide strong individual constraints on 9^ and (3 
(Grego et al. 2001), so we fix {3 at the median value of 0.7 for our sample. Values of 13 in our 
sample are almost entirely in the range 0.5 to 1.0, with most clustered between 0.6 and 0.8; 
this range is also consistent with many independent studies using other cluster samples (e.g., 
Jones & Forman 1999; Mohr et al. 1999; Ettori et al. 2004). To test the effect of changing 
j3, we have repeated the analysis with (3 fixed at 0.6 and 0.8, finding that all changes in the 
parameters are small relative to the 68% statistical uncertainties. We refer to this model as 
the SZE-only model in the remainder of this work. 



4.2. Model Fitting and Likelihood Analysis 

We determine the best-fit values and confidence intervals of the model parameters using 
likelihood analysis based on a Markov chain Monte Carlo (MCMC) method. We work with 
log likelihoods for both the spatial and spectral data. For spatial fitting to X-ray and SZE 
images, the log likelihoods can be written (Cash 1979; Reese et al. 2002) 

\n{Cx-ray) = Yl [A M^i) -Mi- ln(AO] > for X-ray data (Poisson), (9) 

i 

HCsze) = 5^ -- {^^1 + ^^f) Wi for SZE data (Gaussian), (10) 



where Mi and Di are the model prediction and data in pixel i in the X-ray data, while Ai?^ 
and Alf are the difference between the real and imaginary components of the SZE data and 
model at each point i in the Fourier plane and Wi = l/crf is a measure of the Gaussian noise. 
Calculation of the spectral log likelihood is done with XSPEC. The spatial and spectral log 
likelihoods are then added together to compute the joint log likelihood. Bonamente et al. 
(2006) describe the likelihood analysis for the non-isothermal double /9-model; in that case 
the X-ray and SZE central densities were linked and D^ was a free parameter in the fit. Our 
approach differs only in that D^ is fixed for an assumed cosmology and the X-ray and SZE 
central densities are fit separately so that the individual mass results can be compared. We 
therefore describe the likelihood calculation here only in the context of the 100 kpc-cut and 
SZE-only models. 

In the 100 kpc-cut case, the X-ray model is given by Eq. [5] plus a constant X-ray 
background, B^. The model has five free parameters, including the two-dimensional cluster 
position {xc, Vc), central surface brightness S^o, and the shape parameters 6c and (3. The 
background, B^, is determined from the data using a method described in Bonamente et al. 
(2004) and held constant during the model fitting. The model is calculated at each pixel 
position and then multiplied by the exposure map before being compared to the image. X- 
ray point sources and any obvious substructure are masked out, and the masked pixels are 
ignored in the likelihood calculation. Any pixels within 100 kpc of the cluster center are also 
ignored. The Poisson log likelihood is then computed using Eq. [9]. 

The interferometric SZE data are analyzed directly in the Fourier plane, where the 
noise characteristics and spatial filtering of the interferometer are well understood. In both 
the 100 kpc-cut and SZE-only cases, the composite SZE model consists of the isothermal 
/9-model and any point sources detected in the field. The parameters of interest are the 
two-dimensional cluster position, central decrement ATq, 6c, and f3 and the positions and 
fiuxes of the point sources. The SZE /9-model plus point source model are constructed in the 
image plane, multiplied by the primary beam, and fast Fourier transformed to the u-v plane 
for comparison with the data. The Gaussian log likelihood is then calculated using Eq. [10]. 

X-ray spectral likelihoods are calculated using x^ information from XSPEC. A photo- 
absorbed Raymond-Smith thermal plasma model (Raymond & Smith 1977) is fit to the 
0.7 — 7.0 keV X-ray spectral data, with HI column density fixed at the value from Dickey 
& Lockman (1990) and solar abundances from Feldman (1992). XSPEC is used to create a 
table of Tc and metallicity relative to solar, Z, values versus log likelihood for the spectral 
model using the relation \nC = —\x^- 

Best fit values and confidence intervals for model parameters are determined using a 
MCMC method which efficiently handles the large number of parameters involved in the 
model fitting. Implementation of a MCMC method for determining the angular diameter 
distance to the cluster A611 from SZE and X-ray data is described in detail in Bonamente 
et al. (2004), including choice of parameter support and convergence and mixing tests. 
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For the 100 kpc-cut model, we use the hkehhood analysis based on the MCMC method 
and fit the SZE and X-ray data jointly; 9,. and {3 are linked between X-ray and SZE datasets 
(Reese et al. 2000), and the X-ray spectral model is also included. Since the SZE data do 
not constrain 6c and /3 well individually, the X-ray dataset drives the fit to these parameters 
while ATq is allowed to find its best-fit value. The central 100 kpc are excluded from the 
X-ray data but not the SZE data. Since the datasets are independent, the X-ray and SZE 
spatial log likelihoods and spectral log likelihood are added to determine the joint likelihood 
for this model. 

For the SZE-only model, the MCMC method is used to fit the SZE data to an isothermal 
/9-model; (3 is fixed at 0.7 while 9^ and ATq are allowed to find their best fit values. An X-ray 
spectral model with central 100 kpc excluded from the X-ray data is also included. Since the 
SZE and X-ray datasets are independent, the SZE spatial log likelihood and X-ray spectral 
log likelihood are added to determine the joint likelihood for this model. 



4-2.1. Gas Mass, Total Mass and Gas Mass Fraction 

With the best-fit ICM model and X-ray temperature in hand, it is straightforward to 
compute the gas mass and total mass of the cluster. For the /3-model, the enclosed gas mass 
is obtained by integrating the best-fit 3D gas density profile: 

rr/D^ / n2\ -3/3/2 

where A = Annjie^mp Dj^, and /Xg, the mean molecular weight of the electrons, is determined 
from the X-ray spectral data. Assuming the isothermal gas temperature, we can compute 
the gas mass independently from both the X-ray and SZE datasets. For the X-ray data, the 
model central electron density Ue^ can be expressed analytically as (Birkinshaw et al. 1991): 

For the SZE, the model central electron density can be expressed as (e.g., Grego et al. 2001): 

'" \f(.,T.)ToM,T-rkBT,D,,W^/-'rCll3-\)eJ ' 

For the 100 kpc-cut model, we compute the gas mass using Eq. [11] by extrapolating the 
model into the cluster centers. For the non-isothermal /3-model, the gas mass is obtained 
by integrating the best fit central density over a distribution similar to Eq. [11], but which 
accounts for the additional density component. In this case /Xg is treated as a constant, as 
its value changes by only 0.3% over the radial range considered. In the SZE-only analysis, 
the gas mass is computed using Eq. [11] with model central density from Eq. [13]. 
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The total mass, Mtotai, can be obtained by solving the hydrostatic equilibrium equation 

as: 

Mtotai r = -7^ ^e r — ^ + ne r — ^ . 14 

Gfitotrrip \ dr dr J 

Under the isothermal assumption, this reduces to the simple analytic form (e.g., Grego et al. 
2001): 

Mtotai(r) = ^^^, (15) 

which can be used to calculate total masses for both the 100 kpc-cut and SZE-only models. 
For the non-isothermal /3-model, the temperature derivative of Eq. [14] is simple to compute 
numerically. We then compute X-ray and SZE gas mass fractions as /gas = Mgas/Mtotai for 
the sample of 38 clusters. 



4.3. Comparison of the Density Models Fit to X-ray Data 

We now compare the results of the X-ray surface brightness modeling and the cluster 
parameters extracted using different ICM models. The primary goal of this comparison is to 
assess the effects of the isothermal assumption and different treatments of the cluster core. 

Figure 1 shows the isothermal /3-model, non-isothermal /3-model, and 100 kpc-cut model 
as fit to both the non cool-core cluster A1995 and the cool-core cluster A1835. The X-ray 
surface brightness profiles are background subtracted and constructed using concentric an- 
nuli centered on the cluster (e.g., Bonamente et al. 2006). The model fitting is done to the 
entire two-dimensional cluster image. Table 2 lists the spatial {dc and /3), spectral (isother- 
mal spectroscopic gas temperature Tx), and inferred (Mgas, Mtotai, /gas) cluster quantities 
obtained using these three ICM models. The radius r25oo is computed using parameters from 
the 100 kpc-cut models. Temperature profiles and data points for the non-isothermal models 
can be found in Bonamente et al. (2006), who demonstrate that the spectroscopic data are 
well fit by the model temperature profiles. 

In the case of A1995 (a non cool-core cluster), we find that the results are largely in- 
sensitive to the chosen model. The surface brightness profiles appear well fit by all three 
models, and the derived gas masses, total masses, and gas mass fractions are in good agree- 
ment. Although 9c and 13 differ slightly, these parameters are sufficiently degenerate that 
the difference has a very small effect on the derived masses. These results illustrate that 
the simple 100 kpc-cut model works as well as the more sophisticated non-isothermal double 
/3-model for the non cool-core clusters. 

Analysis of the cool-core cluster A1835, on the other hand, highlights the importance of 
a proper treatment of the cluster core. The most striking differences are the masses derived 
from the isothermal /5-model versus those from the other two models. This illustrates the 
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shortcomings of a brute force application of the isothermal /5-model to cool-core clusters. The 
mass discrepancies can be attributed to a poorer fit to the surface brightness at r > 100 kpc; 
this arises because the fit is driven mainly by the extremely high signal to noise data in 
the cluster core. A1835 and other cool-core clusters tend to have extremely small core radii 
and 13 parameters when fit by a single isothermal /3-model. We also find that the X-ray 
spectroscopic temperature, Tx, is biased low when the core is included in the determination 
of Tx, which has an additional impact on the total mass estimate of the cluster. 

The models can be quantitatively compared using a goodness of fit analysis. Goodness 
of fit for the X-ray data is assessed using Monte-Carlo simulations, following Winkler et al. 
(1995). For a given cluster, we construct the best-fit model and compare with the data to 
determine the fit statistic from Eq. [9]. Poisson noise is then randomly added to the best-fit 
model, creating 10^ new realizations, and the fit statistic is calculated for each by comparing 
them with the original best-fit model. The fraction of simulations that give a lower fit 
statistic than that of the best-fit model compared to the data is called the "goodness" (e.g., 
Jonker et al. 2005), with values near 0.5 indicating a good fit to the data and values near 
zero or one indicating a poor fit. 

In the case of A1995, all three models (standard isothermal /9, 100 kpc-cut, and non- 
isothermal) provide acceptable descriptions of the data, with respective goodness values of 
0.416, 0.427, and 0.496. Acceptable goodness values of 0.547, 0.548, and 0.603, for the 
same three models, are also found for A1835. However, the goodness value for the standard 
isothermal /9-model may be biased by the extremely high count rate in the core of this 
cluster. To investigate this, we use the best-fit parameters from the standard /9-model, but 
test them with the 100 kpc-cut dataset. The result is a goodness of zero; the standard (3- 
model is a poor description of the cluster gas distribution between 100 kpc and r25oo- Tests 
on clusters in our sample indicate that 90 — 95% of the cluster gas mass within r25oo lies 
beyond a 100 kpc radius (LaRoque 2005). It is therefore critical that the model give a good 
description of the data between 100 kpc and r25005 which the standard isothermal /3-model 
fails to do. By contrast, when the best-fit non-isothermal /5-model parameters are fit to the 
data between 100 kpc and r25oo, the goodness is 0.776, implying a much better description 
of the data in this region. 



4.4. Testing the Model with Simulations 

The ability of the 100 kpc-cut model to recover the gas mass can be further assessed by 
fitting the model to projected data from simulated clusters with precisely known gas masses. 
This test is carried out using cosmological simulations of galaxy clusters generated with 
the Adaptive Refinement Tree (ART) N-body-|-gasdynamics code (Kravtsov 1999; Kravtsov 
et al. 2002). These simulations include radiative cooling and UV heating, star formation, 
and stellar feedback processes in addition to the standard gas dynamics. Mock Chandra 
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observations of these simulated clusters are performed by using the true gas density, temper- 
ature and metallicity in each gas cell along with the MEKAL spectral emissivity code (Mewe 
et al. 1985) to generate an energy-dependent X-ray flux between 0.1 and 10 keV. This flux 
is then projected along the line of sight, converted to a photon count rate appropriate for 
z = 0.01, and convolved with an instrumental response simulating that of ACIS-I between 
0.7 and 7.0 keV. This generates a counts image similar to the Chandra level 2 event files 
described in §3. The spectroscopic temperature over a region of interest is determined using 
an absorbed MEKAL model in XSPEC (Nagai et al., in preparation). Temperatures derived 
using the MEKAL model differ by at most a few percent from those derived from the rel- 
ativistic Raymond-Smith model used in our analysis. These temperatures should therefore 
well approximate the measured Chandra temperatures for the clusters in our sample. 

We have fit /9-models to the counts images described above for four different simulated 
clusters, spanning a range in mass from 8.7 x 10^^ to 6.0 x 10^^ Mq at r25oo and including 
both relaxed and unrelaxed clusters. Three fits are performed for each cluster; one each 
to projections along three perpendicular axes. The likelihood analysis follows the methods 
of §4.2, and the central 100 kpc of each simulated cluster are excluded from the fit. Any 
obvious gas clumps are also removed from the mock images prior to fitting. A spectroscopic 
temperature is measured between 100 kpc and r25oo! and the masses are determined at r25oo- 
Goodness values for 100 kpc-cut fits to the simulated data are consistently in the range 0.35 
to 0.65, implying an acceptable description of the data. 

Figure 2 shows the mass results for both a highly relaxed and a highly unrelaxed cluster. 
In both relaxed and unrelaxed clusters, the gas mass is recovered to better than 10% between 
0.5r25oo and 2r25oo5 and to better than 5% within r25oo for the twelve different fits. This 
shows that the simple isothermal /9-model with central 100 kpc excluded does an excellent 
job of recovering the true three-dimensional gas mass distribution from the two-dimensional 
projected data at radii encompassing r25oo in simulated clusters. 



5. Additional Sources of Uncertainty 

Before presenting the results, we discuss possible sources of both statistical and sys- 
tematic uncertainty in our measurements and estimate their effects on the X-ray and SZE 
derived gas mass fraction in clusters. Bonamente et al. (2006) address sources of statistical 
and systematic uncertainty in the X-ray and SZE measurements as they apply to calculation 
of the angular diameter distance Da- Some are found to have a negligible effect on our mea- 
surements, such as uncertainties in the assumed Galactic N^ used to calculate spectroscopic 
temperatures, and interference associated with other anisotropics in the CMB. However, we 
do review non-negligible effects discussed in Bonamente et al. (2006) in the context of the 
/gas calculations, and we address additional uncertainties unique to the /gas analysis. The 
effect of these uncertainties on the gas mass fractions is summarized in Table 3. The effects 
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are quantified below in terms of tlieir impact on tfie gas mass fraction for individual clusters. 
Note, however, that in mean results for the sample presented in §6.1, statistical uncertain- 
ties average down by a factor of one over square root of the sample size, while systematic 
uncertainties do not average down. 



5.1. Selection Effects 

Clusters chosen for SZE observation are generally selected on the basis of X-ray luminos- 
ity, Lx- Selecting by X-ray brightness could bias the SZE sample toward elliptical clusters 
elongated along the line of sight; selecting instead by luminosity alleviates this problem. 
However, selecting by L^ may introduce other biases because of the preferential selection 
of cool-core clusters or clusters undergoing major merging events. Another SZE selection 
bias involves point sources; preferentially selecting clusters without bright (> 20 mJy) point 
sources in the field may bias our sample toward non cool-core clusters, or toward clusters at 
higher redshifts whose properties might differ from nearby clusters (e.g., Grego et al. 2001; 
Reese et al. 2002). X-ray selection biases could also be a factor, as we only include clusters 
which can be fit with the 100 kpc-cut and non- isothermal double /3-models. Three clusters, 
CL J0152.7-1357, Abell 520, and MACS J0717.5+3745, are not included in our sample due 
to their highly irregular morphologies. 

We assess the bias associated with sample selection by attempting to construct a sub- 
sample largely free of these effects. The subsample consists of clusters from the BCS catalog 
(Ebeling et al. 1998), which is compiled with data from a single instrument {ROSAT) and has 
a well-defined flux limit of 4.4 x 10"^^ g^g cm'^ s"^ in the R0SATQ.I—2A keV band. The 
subsample contains only clusters with redshifts z < 0.3 to minimize biases associated with 
cluster redshift evolution. To reduce selection effects associated with line-of-sight elongation, 
we impose a flux cut of 6.6 x 10^^^ erg cm^^ s^^ on the subsample, which is 50% higher 
than the BCS flux limit. We further select clusters with intrinsic luminosities in the range 
0.35 < L^ < 1.5 X 10^^ erg s~^ in the 2 — 10 keV Chandra band. This eliminates clusters in 
the extreme high and low luminosity ranges of our sample, which may have unusual dynam- 
ical properties. The final subsample consists of the ten clusters A267, A665, A773, Zw3146, 
A1413, A1914, A2204, A2218, A2261, and R2129 (A1413 and R2129 are not included when 
the SZE-only model is considered, for reasons discussed in §6.1). 

To quantify the effects of selection biases, we calculate the weighted mean /gas values 
using the 100 kpc-cut, non- isothermal, and SZE-only model results (see §6.1) for both the 
full cluster sample and the unbiased subsample. We find that the mean /gas changes by less 
than 3% between the full and unbiased samples, for all three of these models. We therefore 
conclude that selection effects do not appear to have a substantial impact on the gas mass 
fraction or cosmological results. 
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5.2. Kinetic SZE 

The kinetic SZE is a distortion in the CMB spectrum caused by the pecuhar velocity of a 
cluster along the line of sight. Reese et al. (2002) report that, for a cluster with T^ = 8.0 keV 
and with typical line of sight peculiar velocity of 300 km s~^ (Watkins 1997; Colberg et al. 
2000), the kinetic SZE would be 4% of the thermal SZE for our 30 GHz observations. This 
effect can be positive or negative depending on whether the cluster is moving toward us or 
away from us. Since the SZE-derived /gas is linearly proportional to the SZE decrement, the 
kinetic SZE introduces a ±4% additional statistical uncertainty into the /gas measurement 
of each cluster. 



5.3. Radio Point Source Contamination 

Point sources unaccounted for in the SZE data can cause the decrement to be "filled 
in" by the point source flux, diminishing the magnitude of the SZE signal. It is also possible 
to overestimate the magnitude of the decrement if the point source is in a negative beam 
sidelobe. Point sources are identifled in the data using images from DIFMAP. All cluster 
flelds are also cross-referenced with 1.4 GHz maps from the NRAO VLA Sky Survey (NVSS) 
(Condon et al. 1998), which occasionally allow us to distinguish faint point sources with flat 
or gently declining spectra from the noise in the 30 GHz maps. For cases in which NVSS 
shows a bright source that does not show up at 30 GHz, we place a point source model at 
the NVSS position and marginalize over its flux. Strongly inverted sources with 30 GHz 
fluxes near the noise limit would not show up in NVSS, but should be rare. Cooray et al. 
(1998) report spectral indices for 55 point sources with fluxes measured between 1.4 GHz 
and 30 GHz, flnding only three with mildly inverted spectra. 

NVSS contains only point sources with fluxes greater than 2 mJy at 1.4 GHz. To test 
the effect of point sources too faint to be cross-correlated with NVSS, LaRoque et al. (2002) 
randomly distribute point sources in the fleld of the cluster A2163 according to a num- 
ber count versus flux density relationship calibrated from 41 clusters observed with BIMA. 
Considering undetected point sources with intrinsic fluxes between 10 /iJy and 300 fiJj at 
30 GHz, they flnd that the best-flt SZE decrement can change by ±4%. 

Most clusters have massive central cD galaxies, and many of them are radio bright. 
LaRoque et al. (2002) include a point source model at the optical position of the cD galaxy 
in A2163 and marginalize over its flux. They flnd an increase in magnitude of the best-flt 
central decrement of 2%, and the best-flt point source flux is consistent with zero. Reese 
et al. (2002) uses a similar treatment, placing point source models at the centers of the 
clusters A697, A2261, and A1413. The best-flt fluxes are all consistent with zero, and again 
increase the magnitude of ATq by only ~ 2%. We conclude that the effect of undetected 
point sources is likely dominated by the off-center sources, and should add an additional 
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±4% statistical uncertainty to the SZE gas mass fractions for each cluster. 



5.4. Instrument Calibration 

The absolute calibration of the ACTS response is uncertain at the 5% level in the 0.7 — 
7.0 keV band (http://cxc.harvard.edu/cal), after CTI corrections and correction for con- 
tamination on the optical blocking filters are applied. Measurements of Tx are also subject 
to ~ 5% systematic uncertainty, owing mainly to calibration errors between energy bins. 
The Chandra calibration uncertainty therefore translates to ~ 3% uncertainty in the X-ray 
gas mass {M^^ ^ oc S^o^''^) and 5% in the HSE total mass (Mtotai oc Tx) for a 6% uncer- 
tainty in the X-ray /gas- Absolute calibration of the BIMA/OVRO interferometric imaging 
system is known conservatively to 4% (Reese et al. 2002), which when combined with the 
uncertainty in Tx translates to a 6% systematic uncertainty in the SZE derived gas mass 
(M oc ATq/Tx) and 8% systematic uncertainty in the SZE-derived /gas- 



5.5. X-ray Background 

X-ray background is usually estimated from regions as far from a cluster as possible; 
however, due to space constraints on the chip, the background regions might contain some 
amount of cluster flux, particularly for low- 2; clusters which subtend a larger solid angle on 
the sky. We therefore test the effect of overestimating the X-ray background due to cluster 
contamination on a relatively nearby cool core cluster, Zw3146. This cluster is on chip 
13, so the original background is extracted from circular regions on chips 10, II, and 12 at 
appropriate distances from the readouts (e.g., Bonamente et al. 2004). To assess the effect 
of overestimating the background, we move the extraction regions closer to the cluster until 
the best-fit background increases by ~ 25% due to cluster contamination. This should be 
seen as an extreme case, only applicable to a small number of bright, nearby clusters in our 
sample. The spectral fits are redone with these new background regions, and the resulting 
temperature and higher background level are applied in the X-ray/SZE joint fit. 

We find that the higher X-ray background modifies the best-fit parameters, but the 
magnitude of the effect on individual parameters is small (< 5%) with a combined effect 
on /gas of ~ —7%. We repeat this same analysis on one of the non cool-core clusters, 
A2259, which has a much smaller central surface brightness to background ratio, and find 
that the results are much the same. Given that typical statistical uncertainty in B^ is at 
most ~ 5%, much less than the systematic uncertainty of 25% considered here, fixing the 
background in our analysis should have a negligible effect on the /gas results presented in 
this paper. Accounting for the fact that some clusters are more susceptible than others to 
an overestimated X-ray background, we add a one-sided +2% systematic uncertainty to the 
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X-ray derived /gas results 



5.6. Effects of Asphericity 

All cluster analysis presented in this paper uses a spherical model. Here we assess the 
effects of this assumption on the results presented in this paper. To assess this effect, we refit 
the X-ray data of all 38 clusters using a 100 kpc-cut isothermal /9-model and accounting for 
projected elhpticity rj and position angle of the clusters on the sky. We find a mean projected 
ellipticity (77) = 0.79 with rms scatter 0.09, where rj is the ratio of a cluster's minor axis 
to major axis. These results are consistent with the previous studies of Mohr et al. (1995). 
Accounting for the inclination angle (e.g., Cooray 1998), an ellipsoidal cluster with rj = 0.8 
and axis of symmetry inclined at 45° to the line of sight can have a derived total mass in 
error by 12% and derived gas mass in error by 6% if modeled as a sphere. Extending this to 
the more complex triaxial distributions, we conservatively estimate a statistical uncertainty 
of 10% in the X-ray derived Mgas and 15% in Mtotai- This translates to a ~ 20% statistical 
uncertainty in the X-ray derived gas mass fractions, given that the changes in gas mass and 
total mass can either add or partially cancel depending on inclination angle and geometry 
(Piffaretti et al. 2003). The effect on the SZE /gas results should be smaller since the SZE 
depends linearly on density. Grego et al. (2000) calculate the gas mass fraction of A370 
assuming both prolate and oblate geometries and find a resulting uncertainty in /gas of 20%. 
This cluster has rj = 0.64, however, so we assume a milder 10% statistical uncertainty in the 
SZE /gas results for each cluster. 



5.7. Hydrostatic Equilibrium 

The HSE condition is a key assumption that enables measurements of the gravitationally 
bound mass of clusters from the X-ray and SZE data. The assumption is that the ICM is 
supported against gravitational collapse only by thermal pressure of the hot intracluster gas. 
The hydrostatic mass would therefore underestimate the true gravitational bound mass of 
the cluster if non-thermal pressure support is present. One form of non-thermal pressure 
comes from turbulent gas motions in clusters. Recent numerical simulations suggest that 
this provides about 10% of the total pressure support in clusters (Faltenbacher et al. 2005; 
Rasia et al. 2006; Lau et al. 2006) and similar results are obtained using XMM-Newton 
observations of the Coma cluster (Schuecker et al. 2004). More direct measurements of the 
turbulent gas motions in clusters, however, await measurements of Doppler broadening of 
the iron lines with high-resolution X-ray spectroscopy experiments (Inogamov & Sunyaev 
2003). Cosmic rays and magnetic fields may also provide non-thermal pressure support. 
While direct measurements of the non-thermal pressure support may be difficult, detailed 
comparisons of the hydrostatic mass with other mass estimates for the same cluster (e.g.. 



gravitational lensing and the velocity dispersion of galaxies in clusters) will likely provide 
an important check on the hydrostatic mass estimate and potentially interesting constraints 
on the non-thermal pressure support in clusters. For now we assign a one-sided systematic 
uncertainty of —10% to the total masses which accounts for the assumed contribution from 
non-thermal pressure. We caution, however, that this estimate is somewhat uncertain and 
this effect will need to be investigated in greater detail as more data become available. 



5.8. Isothermal Assumption 

The isothermal assumption can potentially affect the gas mass fraction measurements 
and the cosmological constraints through its effects on both gas mass and total mass. The 
X-ray derived gas mass is affected little by the isothermal assumption, because the X-ray 
emission has a very weak dependence on the assumed temperature (A oc T^/^). The SZE 
derived gas mass is more sensitive to this assumption since it depends linearly on the cluster 
gas temperature. Our choice of a radial range between 100 kpc and r25oo should minimize the 
bias due to the temperature structure for two reasons: (1) use of the X-ray temperature with 
the 100 kpc core exclusion should avoid the most significant bias associated with cool-core 
clusters, and (2) focusing on the inner r25oo should minimize the bias due to the temperature 
decline observed in the outskirts of clusters (Vikhlinin 2005). 

The total mass derived from the X-ray and SZE data is also affected by the isothermal 
assumption. If the temperature profile deviates from the isothermal gas distribution, the 
average temperature depends on the weighting scheme. For example, the spectroscopic 
temperature, T^pec, may differ by about 10% from the gas mass-weighted temperature T^g 
in the radial range between 70 kpc and rsoo (Vikhlinin 2005). But the difference is likely 
smaller in the radial range (100 kpc < r < r25oo) considered in this study. We therefore 
expect that the bias on the cluster masses and hence our gas mass fraction measurements is 
not significant. 

The effect of the isothermal assumption can be quantified using the results of Table 7, 
presented in §6.1 below. There we find that the X-ray and SZE weighted mean gas mass 
fractions are nearly identical for the non- isothermal model, while the SZE gas mass fraction 
is ~ 5% larger than the X-ray value for the 100 kpc-cut model. This hints at a possible 
systematic bias of about —5% in the spectral temperature under the isothermal assumption, 
since the SZE /gas has an additional T^^ dependence relative to the X-ray. We therefore 
assign an additional —5% systematic uncertainty to the gas mass fractions derived from the 
X-ray 100 kpc-cut fits, and a —10% systematic uncertainty to the gas mass fractions from 
the SZE 100 kpc-cut fits. 
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5.9. Cluster Substructure 

The X-ray derived gas mass fractions are particularly sensitive to any non-negligible 
small scale structure in the intracluster gas. When gas clumps, it has density higher than 
the local average, (ng), and the X-ray surface brightness is increased by a factor of C = 
(rig) / (rie) relative to the SZE decrement (Nagai et al. 2000; Mathiesen & Evrard 2001). 
Since Ug oc Sj for the X-ray data and n^ oc AT for the SZE data. X-ray gas mass fractions 
would be systematically overestimated by a factor of C^^"^ with respect to SZE gas mass 
fractions in the presence of significant clumping. Cluster substructure has two important 
effects on the X-ray observables. First, gas clumps can bias the measured X-ray gas mass 
high by C^/^. Second, cluster substructure can bias the average spectroscopic temperature 
of the cluster as a whole (e.g., Mathiesen & Evrard 2001; Mazzotta et al. 2004). 

The Chandra observations provide important handles on both of these effects, because 
their superb resolution and sensitivity enable detection of the prominent clumps that con- 
tribute significantly to the X-ray surface brightness and the determination of the average 
temperature of the cluster. We test this using the mock Chandra observations of the sim- 
ulated clusters introduced in §4.4. The test shows that the gas mass is recovered to better 
than 5% using the simple 100 kpc-cut model, indicating that the clumping factor should be 
small if the most prominent clump is removed from the image. These results should also 
hold for the high z clusters, as a significant fraction of the low energy photons from the cool 
clumps are redshifted out of the 0.7 — 7 keV band, making the instrument less sensitive to 
clumping effects (Nagai et al., in preparation). 

For the same reasons, the bias on the Chandra spectroscopic temperature is also small. 
Tests on several clusters in our sample show that not removing bright clumps may bias the 
temperatures low by 10 — 15%, in reasonable agreement with the previous findings (Mathiesen 
& Evrard 2001; Mazzotta et al. 2004). Further tests, leaving in only fainter clumps, show a 
decrease in spectroscopic temperature of only 2 — 3%, smaller than the assumed bias due to 
the isothermal assumption. 

We therefore conclude that there is no evidence that cluster substructure significantly 
affects our results when bright spots are removed from the X-ray images prior to fitting. 
The agreement in the mean X-ray and SZE gas mass fractions, particularly for the non- 
isothermal model which avoids the bias due to the isothermal assumption, is consistent with 
this conclusion (cf. §6.3). 
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6. Results and Discussion 

6.1. Gas Mass Fraction Results 

In Table 4 we present the X-ray and SZE measurements of cluster gas mass fractions as 
determined from the 100 kpc-cut model (§4.1.1). The X-ray and SZE data are fit jointly, with 
the cluster core radius and 13 parameter linked between datasets while the central surface 
brightness Sxo and central SZE decrement ATq are allowed to vary. The gas mass fractions 
are plotted as a function of redshift in Figure 3. X-ray surface brightness profiles and the 
best-fit 100 kpc-cut models are shown in Appendix A for the 38 clusters in our sample. 
The X-ray surface brightness profiles are well described by the isothermal /9-model beyond 
100 kpc from the cluster center. 

In Table 5 we present X-ray and SZE gas mass fractions determined from the non- 
isothermal model (§4.1.2). The X-ray and SZE data are again fit jointly with shape pa- 
rameters linked, while the X-ray and SZE central density parameters are allowed to vary 
separately. The gas mass fractions are plotted as a function of redshift in Figure 3. Surface 
brightness profiles for the non-isothermal model are shown in Bonamente et al. (2006), and 
show consistently good fits to the cluster radial profiles out to r25oo- 

In Table 6 we present SZE gas mass fractions determined from the SZE-only model 
(§4.1.3). We find that this model constrains ATq and 9^ individually in 23 of our 38 clusters, 
and those 23 are shown in the table. Gas mass fractions derived with this method are plotted 
versus redshift in Figure 3. 

The mean X-ray and SZE gas mass fractions for all three models are summarized in 
Table 7, where the uncertainties are statistical followed by systematic at 68% confidence. The 
statistical uncertainties include the effects of cluster asphericity, radio point sources, and the 
kinetic SZE, combined in quadrature; systematic uncertainties include errors associated with 
instrument calibration. X-ray background, and the assumptions of HSE and isothermality 
(see §5 and Table 3). The results show that the X-ray and SZE measurements of the mean 
gas mass fraction in clusters agree to better than 10% for the three ICM models considered. 



6.2. Effects of Cluster Cores on Derived Gas Mass Fractions 

The results summarized in Table 7 show that the X-ray and SZE derived gas mass 
fractions agree remarkably well as long as care is taken to account for the effect of the 
cluster core on the X-ray emission. The results show that the core can be satisfactorily 
accounted for by either ignoring the core in fits to the X-ray data (the 100 kpc-cut model) 
or modeling the ICM with the non-isothermal double-/? model. 

It is also noteworthy that the gas mass fraction results derived from spatial fits to only 
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the SZE data including the core also agree remarkably well with the above models. This 
indicates that the SZE is not very sensitive to the cluster core, and is in line with expectations 
since the SZE is a probe of the ICM pressure. 

The effect of cluster cores on the gas mass fractions is further investigated by comparing 
results for cool-core and non cool-core subsamples. There is a population of clusters in our 
sample with bright and sharply peaked cores (cf. Appendix A), including A586, MC0744, 
Zw3146, A1413, MC1311, A1689, R1347, MS1358, A1835, MC1423, A2204, A2261, and 
R2129. These clusters have central cooling times (calculated with the cores included in the 
X-ray data) tcooi ~ skBTxnH jggg ^j-^g^]^ O.StHubbie, and we refer to this subsample as the 
cool-core sample because the sharply peaked X-ray emission is indicative of strong radiative 
cooling in the cluster core. The results of a comparison of cool-core and non cool-core 
weighted mean gas mass fractions are shown in Table 8, where in addition to the three 
main ICM models considered above, we include results from two additional naive models 
which make no effort to compensate for the cluster cores: (1) a no-cut isothermal /9-model 
fit jointly to the full X-ray and SZE datasets, and (2) a spatial model fit to only the SZE 
data as before, but using the X-ray temperature derived from all of the X-ray data instead 
of the data with central 100 kpc excluded. 

Apart from a small < 10% offset in the X-ray results, the derived gas mass fractions for 
the 100 kpc-cut and non-isothermal models for the cool-core and non cool-core clusters are 
in excellent agreement. Figure 3 also demonstrates this agreement, as there is no systematic 
offset seen between the cool-core (triangles) and non cool-core (squares) clusters. The results 
for the naive models clearly show the importance of accounting for the core in the cool- 
core subsample; large offsets are found for the derived gas mass fraction of the cool-core 
subsample, while the same naive models provide perfectly consistent results for the non 
cool-core subsample. The primary cause of the biased results of the naive models for the 
cool-core clusters is the poor spatial fits and the low Tx- A depressed Tx results in a lower 
derived total cluster mass as well as a higher SZE gas mass. 



6.3. Constraints on Clumping in the Intracluster Gas 

The comparison of the X-ray and SZE gas mass fraction measurements shows that the 
results are in good agreement. These results suggest a clumping factor of the ICM consistent 
with zero. The high resolution Chandra data allow one to remove any visible clumps before 
the data are fit, and the agreement between X-ray and SZE gas mass fraction measurements 
indicates that any clumps below the noise level of the data have a negligible effect on the 
mass measurements. Further indications in favor of a negligible clumping factor are discussed 
in 55.9. 
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6.4. Comparison with previous studies 

Finally, we compare our X-ray and SZE gas mass measurements to the results from 
previous studies in the literature. We first compare our X-ray results with Chandra results 
obtained using a different analysis technique (Allen et al. 2002, 2004)^, in which the cluster 
mass is modeled assuming that the cluster gas is in hydrostatic equilibrium with a NFW 
potential and using a deprojected 3D gas temperature profile. These authors report mean 
gas mass fractions of /gas = O.f 13 ± 0.005 for a sample of six highly relaxed clusters (Allen 
et al. 2002) and /gas = 0.117 ± 0.002 for a larger sample of 26 relaxed clusters (Allen et al. 
2004); these results are in good agreement with our own. We also compare the results of 
eight clusters we have in common with the Allen et al. (2004) dataset and find that the 
measurements of individual clusters are consistent with no systematic offset. 

For the SZE gas mass fractions, we compare results of our SZE-only model fitting 
directly with the results of Grego et al. (2001). There are many similarities between our 
cluster sample and analysis method and those of Grego et al.; for instance, they use BIMA 
and OVRO data, and 13 of our 23 clusters have data in common with their sample. They also 
use an isothermal /9-model to analyze the SZE data. Their approach also differs somewhat 
from ours, in that they (1) use spectroscopic temperatures from the literature (generally 
from ASCA)^ and (2) measure the gas mass fraction at 65" from the cluster center, then 
scale the results to rsoo using relations calibrated from numerical simulations (Evrard et al. 
1996; Evrard 1997). Despite these differences, Grego et al. measure a mean gas mass fraction 
of /gas = O.OSI^q;!^?^ /i~\ which is equivalent to 0.116lo!oi6 ^^^ ^ = ^-^i i^ good agreement 
with our SZE derived gas mass fraction (cf. Table 7). Similarly, if we calculate gas mass 
fractions from our SZE-only model at 65" and scale to rsoo using the same relations, we 
obtain a weighted mean /gas of 0.086 ± 0.007 /i~^, in good agreement with the Grego et al. 
results. The cluster gas mass fraction results derived from the implied SZE for nearby 
clusters obtained from correlating known clusters with WMAP data is also consistent with 
our results (Afshordi et al. 2005). 



7. Constraints on Cluster Physics and Cosmology 

7.1. The cluster baryon budget 

Currently, the relation between /gas and the cosmo logical baryon fraction VLb/^m is 
not accurately known (Ettori et al. 2006). However, since recent CMB studies provide 



^Although mean gas mass fraction results for other cluster samples are available in the literature (e.g., 
Evrard 1997; Mohr et al. 1999; Ettori et al. 2003; Sanderson et al. 2003), the detailed comparisons with our 
results are non-trivial as they are often calculated at different radii and assume different cosmologies. 
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precise constraints on cosmological parameters (e.g., Spergel et al. 2006), we can use our /gas 
measurements to obtain constraints on the ratio 



/gas 



(fis/fi 



M 



(16) 



Combining our /gas results with the WMAP values of VIb and VIm {^sh'^ = 0.0223±[]:[j[][!^, 
^M^^ = 0.127;*]q'q53) ^^^ the Hubble Space Telescope measurement oi h = 0.72 ±0.08 (Freed- 
man et al. 2001), we obtain the constraints on r^gas given in Table 9, where uncertainties are 
again statistical followed by systematic at 68% confidence. Since the /gas data are calculated 
assuming ^7^=0. 3, r2A=0.7, h = 0.7, very similar to the WMAP constraints, no significant 
bias is expected in the results of Table 9. These results indicate that the ratio of the gas mass 

fraction in clusters within r25oo to the cosmic baryon fraction, 77—77^ — , is 0.68to;}6 where 
the range includes statistical and systematic uncertainties, in accord with the conclusions 
reached by Afshordi et al. (2005). This ratio provides a benchmark for simulations of galaxy 
clusters. We note that recent numerical simulations that include radiative cooling and star 
formation have baryon budgets in line with our results, although they find that the gas tends 
to overcool in the cluster cores (Kravtsov et al. 2005). 



7.2. Cosmological Parameters 

Cluster gas mass fractions can provide strong constraints on cosmological parameters. 
The normalization of the cluster gas mass fraction depends sensitively on the ratio of Qb/^m 
(e.g.. White et al. 1993; David et al. 1995; Myers et al. 1997; Mohr et al. 1999; Ettori & 
Fabian 1999; Grego et al. 2001; Allen et al. 2002), while the redshift evolution of the gas 
mass fraction depends sensitively on Qm, ^a and even the equation of dark energy, w (Sasaki 
1996; Pen 1997; Allen et al. 2004). Methods involving the normalization, however, depend 
sensitively on a priori knowledge of the cluster baryon budget which, as discussed in § 7.1, is 
poorly known. By contrast, methods using the redshift evolution are powerful in that they 
are independent of the cluster baryon budget and even prior knowledge on other cosmological 
parameters such as Qb and h. 

We therefore consider constraints on cosmological parameters from the redshift evolution 
of the cluster gas mass fraction. For parameter estimation, we follow the procedure described 
in Allen et al. (2004), and fit the data to the model. 



)fid qM\ X (3/2 [X-ray], 1[SZE]) 



pmod|' N ^ AT I -^aI^,"m^"A 



where we assume a ACDM cosmology with Qm=0.3, r2A=0.7, as a fiducial cosmology, indi- 
cated by "fid" . Using this model, we fit an arbitrary normalization A^, and grid the variables 
Hm and Ha between and 1.3 with the WMAP constraint Hm + ^a = l-OOSiaoi? (Spergel 
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et al. 2006). We cannot place useful constraints on Qm and ^a from the X-ray or SZE data 
without a prior on the total density parameter. 

Marginalizing over the WMAP total density prior, we find best-fit parameters of VLm = 
0.26l[j:?^, ^A = 0.74+°:^^ from the 100 kpc-cut X-ray data and VLm = OAOtoH and ^a = 
0.601q'5^ from the isothermal P model analysis. Similarly, the non-isothermal model returns 
Qm = 0.40+|]J^, ^A = 0.60+[J|^ from X-ray and ^m = 0.18+[J:?^, J^a = 0.82+°;^^ from SZE. 
Our results are consistent with independent tests, such as those of WMAP, and inconsis- 
tent with a matter dominated universe (^2^=1, ^a=0) model. In Figure 4 we show the 
/gas data, the best-fit model using the WMAP prior (solid line), the best-fit model without 
the WMAP prior (unconstrained; dashed line), and the model corresponding to two other 
fiducial cosmologies (^2^=1, fiA=0: dotted line; ^7^=0. 3, r2A=0: dash-dotted line). Only 
the 100 kpc-cut results are plotted for simplicity. It is evident that the X-ray data disfavors 
the fiat, matter-dominated model. Interestingly, systematic uncertainty should be negligible 
compared to statistical uncertainty with this method, since most of the systematic uncer- 
tainties affect the normalization and do not introduce trends with redshift. Cosmo logical 
constraints from this method are therefore expected to improve as the cluster sample size is 
increased. 



8. Conclusions 

We have presented X-ray and SZE measurements of gas mass fractions for a sample of 
38 massive galaxy clusters using Chandra X-ray observations and BIMA/OVRO SZE radio 
interferometric observations. Three models for the cluster gas distributions are used in the 
gas mass fraction derivations: (1) an isothermal /9-model fit to the X-ray data at radii beyond 
100 kpc and to all of the SZE data, (2) a non-isothermal double /9-model fit to all of the 
X-ray and SZE data, and (3) an isothermal /3-model fit to the SZE spatial data with f3 fixed 
at 0.7 and a spectral temperature from the 100 kpc-cut X-ray data. 

The cluster gas distributions are well fit by the non-isothermal model out to r25oo- 
We have also demonstrated that the simple 100 kpc-cut isothermal model provides a good 
description of the intracluster medium (ICM) outside of the cluster core for objects with a 
wide range of morphological properties. 

For all three models, the mean gas mass fractions determined from the X-ray and SZE 
data are in excellent agreement. Specifically, with the isothermal /9- model we find a mean 
X-ray gas mass fraction /gas(X-ray)= O.llOtg'ogg tg'o^s and a mean SZE gas mass fraction 
/gas(SZE)= 0.1161q;qo5 -ao26' "where uncertainties are statistical followed by systematic at 
68% confidence. For the non-isothermal double /5-model, the mean gas mass fractions are 
/gas(X-ray)= 0.119+°:°°^+°:°?! and /gas(SZE)= 0.121+°:°05 -ao?'- For the SZE-only model, 

f /'Q7"R'l— n lon+0-009 +0.009 
Jgasy^^i^)— U.i^U_goo9 -0.027- 
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It is noteworthy that the gas mass fractions derived from spatial fits to only the SZE 
data including the core also agree remarkably well with the other models. This indicates 
that the SZE is largely insensitive to structure in the cluster core. The lack of sensitivity to 
the core structure is expected since the SZE is a probe of the gas pressure; this is promising 
for extracting cosmology from upcoming large scale SZE cluster surveys (e.g., see Carlstrom 
et al. 2002). 

To investigate the effect of the cluster cores on the results, we divided the cluster sample 
into two subsamples depending on whether or not the clusters exhibit cool cores. The three 
ICM models above lead to consistent gas mass fraction results for each subsample. However, 
when the subsamples are fitted with naive ICM models that make no effort to account for 
the cluster cores, the resulting gas mass fractions for the cool-core subsample show large 
offsets from the above results, while the results for the non cool-core clusters are consistent 
with the above results. This clearly indicates the importance of accounting for the cores in 
clusters that exhibit cool cores. 

Understanding the baryon budget in clusters is one of the most important aspects of 
cluster physics. Our X-ray and SZE measurements indicate that the ratio of the gas mass 

fraction within r25oo to the cosmic baryon fraction, o ,T) — , is O.GS^gJg where the range 
includes statistical and systematic uncertainties at 68% confidence. This result is in accord 
with the conclusions reached by Afshordi et al. (2005). The measurement provides a bench- 
mark for simulations of galaxy clusters. We note that recent numerical simulations that 
include radiative cooling and star formation have baryon budgets in line with our measure- 
ments, although they find that the gas tends to overcool in the cluster cores (Kravtsov et al. 
2005). 

Our results also have implications for cosmological constraints based on measurements 
of cluster gas mass fractions. Constraints from the normalization of the gas mass fraction 
depend sensitively on a priori knowledge of the cluster baryon budget r^gas, which, as dis- 
cussed above, is still somewhat uncertain. However, a cosmological test based on the redshift 
evolution of the gas mass fraction is independent of the cluster baryon budget and even prior 
knowledge of cosmological parameters such as Qb and h. Marginalizing over the WMAP 
total density prior, we use this method and find best-fit parameters of VLm = 0.26t.Q'1l, 
n^ = 0.74+°:^^ and ^m = 0.40+E5:26> ^A = 0.601^:^^ from the 100 kpc-cut X-ray and SZE 
samples, and ^m = 0.40+[J|^, ^a = 0.60+^|^ and Qm = 0.18+^;?^, ^a = 0.82+°;^^ from the 
non- isothermal X-ray and SZE samples, respectively. Our results are consistent with stan- 
dard ACDM cosmology, and inconsistent with ^2^=1, r2A=0. Interestingly, the systematic 
uncertainty for this method is negligible compared to statistical uncertainty, as most sources 
of systematic uncertainty are not expected to introduce trends with redshift. Cosmological 
constraints using this method are therefore expected to improve rapidly as the cluster sample 
size increases with upcoming X-ray and SZE cluster surveys. 
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Fig. 1. — : Comparison of X-ray surface brightness profiles for the isothermal /3-model, non- 
isothermal /3-model, and 100 kpc-cut model as applied to both the non cool-core cluster 
A1995 (upper panels) and the cool-core cluster A1835 (lower panels). The dotted vertical 
line in the 100 kpc-cut panels denotes the 100 kpc radius. 
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Fig. 2. — : Ability of the 100 kpc-cut model to recover the gas mass at r25oo in two simulated 
clusters, one highly relaxed (top panel) and one highly unrelaxed (bottom panel). The upper 
part of each panel shows the true gas mass profile (dark solid line) with dash-dotted lines 
showing profiles recovered from fitting the model to x, y, and z projections. The lower part 
of each panel shows the fractional deviation of recovered mass from true mass as a function 
of radius. 
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Fig. 3. — : Gas mass fraction versus redshift derived from the X-ray data (left panels) and 
from the SZE data (right panels). Triangles (blue) denote cool-core clusters and squares 
(red) denote non cool-core clusters. Dashed lines show the weighted mean /gas for each 
sample. Error bars are statistical from Tables 4, 5, and 6. 
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Fig. 4. — : 100 kpc-cut derived gas mass fractions versus redshift for X-ray (left) and SZE I 
(right) data, showing best fit models for different cosmologies. The solid lines (black) show 
the best fit cosmologies when the prior Qm + ^a = 1-00 is assumed. The dashed (green) hues 
show the best fit cosmologies when no prior is assumed. The dash-dotted (orange) lines show 
the best fit normalization when an f2j\,/=0.3, r2A=0 cosmology is assumed, and the dotted 
(magenta) line shows the best fit normalization for a cosmology with r2j\/=1.0, r2A=0. 



Table 1. Cluster Data 
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13 


20.2 
19.4 


07 44 52.8 


+39 27 26.7 


8 


17 


07 44 52.4 


+39 27 33.2 


LaRoque ct al. (2003) 


ABELL 611 


0.288 


3194 


S3 


36.1 


08 00 56.6 


+36 03 24.1 


- 


57 


08 00 56.5 


+36 03 22.9 


Struble & Rood (1999) 


ABELL 665 


0.182 


3586 
531 


13 
13 


29.7 
9.0 


08 30 58.1 


+65 50 51.6 


52 


16 


08 30 58.6 


+65 50 49.8 


Struble & Rood (1999) 


ABELL 697 


0.282 


4217 


13 


19.5 


08 42 57.5 


+36 21 56.2 


- 


47 


08 42 57.8 


+36 21 54.5 


Struble & Rood (1999) 


ABELL 773 


0.217 


533 

3588 


13 
13 


11.3 
9.4 


09 17 52.8 


+51 43 38.9 


26 


66 


09 17 53.5 


+51 43 49.8 


Struble & Rood (1999) 


ZW 3146 


0.291 


909 


13 


46.0 


10 23 39.7 


+04 11 09.5 


25 


15 


10 23 37.8 


+04 11 17.8 


Allen et al. (1992) 


MS 1054.5-0321 


0.826 


512 


S3 


89.1 


10 56 59.4 


-03 37 34.2 


- 


43 


10 56 59.1 


-03 37 34.0 


Luppino & Gioia (1995) 


MS 1137.5-1-6625 


0.784 


536 


13 


77.0 


11 40 22.3 


+66 08 16.0 


88 


- 


11 40 23.1 


+66 08 05.3 


Donahue et al. (1999) 


MACS J1149.5-I-2223 


0.544 


1656 
3589 


13 
13 


18.5 
20.0 


11 49 35.5 


+22 24 02.3 


39 


— 


11 49 34.9 


+22 23 54.8 


LaRoque et al. (2003) 


ABELL 1413 


0.142 


1661 
537 


13 
13 


9.7 
9.6 


11 55 18.0 


+23 24 17.0 


28 


- 


11 55 17.7 


+23 24 39.5 


Struble & Rood (1999) 


CL J1226.9+3332 


0.890 


3180 
932 


13 

S3 


31.7 
9.9 


12 26 57.9 


+33 32 47.4 


33 


- 


12 26 58.0 


+33.32 57.9 


Ebeling et al. (2001) 


MACS J1311. 0-0310 


0.490 


3258 


13 


14.9 


13 11 01.7 


-03 10 38.5 


39 


- 


13 11 02.2 


-03 10 45.6 


Allen et al. (2004) 


ABELL 1689 


0.183 


1663 
540 


13 
13 


10.7 
10.3 


13 11 29.5 


-01 20 28.2 


16 


26 


13 11 29.1 


-01 20 29.7 


Struble & Rood (1999) 


RX J1347.5-1145 


0.451 


3592 


13 


57.7 


13 47 30.6 


-11 45 08.6 


22 


3 


13 47 30.6 


-11 45 12.3 


Schindler et al. (1995) 


MS 1358.4-1-6245 


0.327 


516 


S3 


48.1 


35 59 50.6 


+62 31 04.1 


70 


- 


13 59 50.2 


+62 31 07.0 


Gioia & Luppino (1994) 


ABELL 1835 


0.252 


495 
496 


S3 
S3 


19.5 
10.7 


14 01 02.0 


+02 52 41.7 


27 


23 


14 01 01.8 


+02 52 45.6 


Struble & Rood (1999) 


MACS J1423.8+2404 


0.545 


4195 


S3 


115.6 


14 23 47.9 


+24 04 42.6 


35 


- 


14 23 47.7 


+24 04 37.3 


LaRoque et al. (2003) 


ABELL 1914 


0.171 


3593 


13 


18.9 


14 26 00.8 


+37 49 35.7 


24 


- 


14 26 01.3 


37 49 38.6 


Struble & Rood (1999) 



00 



35 













Table 1 — Continued 












z 






Chandra X-ray Data 






Interferometric SZE Data 






Cluster 


ObsID 


Chip 


(ks) 


(hh mm ss) 


(o / // ) 


BIMA (hr) 


OVRO (hr) 


(hh mm ss) 


(o / //) 


z reference 






542 


13 


8.1 
















ABELL 1995 


0.322 


906 


S3 


56.7 


14 52 57.9 


-1-58 02 55.8 


50 


58 


14 52 58.1 


+58 02 57.0 


Patel et al. (2000) 


ABELL 2111 


0.229 


544 


13 


10.3 


15 39 41.0 


-1-34 25 08.8 


36 


- 


15 39 40.2 


+34 25 00.4 


Struble & Rood (1999) 


ABELL 2163 


0.202 


1653 


11 


71.1 


16 15 46.2 


-06 08 51.3 


23 


37 


16 15 43.6 


-06 08 46.6 


Struble & Rood (1999) 


ABELL 2204 


0.152 


499 
6104 


S3 
13 


8.6 
9.6 


16 32 46.9 


-1-05 34 31.9 


30 


— 


16 32 46.6 


+05 34 20.6 


Struble & Rood (1999) 


ABELL 2218 


0.176 


1666 


10 


41.7 


16 35 51.9 


+66 12 34.5 


32 


70 


16 35 48.7 


+66 12 28.1 


Struble & Rood (1999) 


RX J1716.4+6708 


0.813 


548 


13 


51.7 


17 16 48.8 


-1-67 08 25.3 


37 


- 


17 16 51.2 


+67 07 49.6 


Henry et al. (1997) 


ABELL 2259 


0.164 


3245 


13 


10.0 


17 20 08.5 


4-27 40 11.0 


25 


- 


17 20 09.0 


+27 40 09.4 


Struble & Rood (1999) 


ABELL 2261 


0.224 


550 


13 


9.1 


17 22 27.1 


-1-32 07 57.8 


23 


40 


17 22 26.9 


+32 07 59.9 


Struble & Rood (1999) 


MS 2053.7-0449 


0.583 


551 
1667 


13 
13 


44.3 
44.5 


20 56 21.2 


-04 37 47.8 


— 


154 


20 56 21.0 


-04 37 47.2 


Stocke et al. (1991) 


MACS J2129.4-0741 


0.570 


3199 
3595 


13 
13 


8.5 
18.4 


21 29 26.0 


-07 41 28.7 


- 


24 


21 29 24.9 


-07 41 43.9 


LaRoquc et al. (2003) 
Ebeling et al. (1998) 


RX J2129.7+0005 


0.235 


552 


13 


10.0 


21 29 39.9 


-(-00 05 19.8 


47 


- 


21 29 38.1 


+00 05 12.4 


MACS J2214.9-1359 


0.483 


3259 
5011 


13 
13 


19.5 
16.1 


22 14 57.3 


-14 00 12.3 


41 


11 


22 14 58.4 


-14 00 10.9 


Note C') 


MACS J2228.5-I-2036 


0.412 


3285 


13 


19.9 


22 28 33.0 


-1-20 37 14.4 


39 


- 


22 28 33.1 


+20 37 14.2 


Bohringer et al. (2000) 



^Redshift derived from the Fe lines in the Chandra x-ray spectrum, this paper. 
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Table 2. Comparison of Parameters for Different ICM Models 



Cluster 


Model 


Oc 


/3 


Tx 


^gas(^250o) 


M;otal(^250o) 


/gas 






n 




(keV) 


(lOi^Mo) 


(lO^^Mo) 




A1995 


Isoth. p 


50.3t}| 


0.92i;n^^ 


8.63+°i^ 


q p:9+0.12 
'-'•'-'^-0.12 


4 Rq+0.39 


0.075+°:°°! 




Non-Isoth. p 


50.8111 


0.925tH^^ 




q 77+0.21 
"J- ' '-0.23 


4 Q7+0.75 
^■<J'-0.71 


0.086t°:°^^ 




100 kpc-cut 


57 5+2.7 


1 01+°-°'' 


8.22+°:^^ 


q e;i+0.14 


A 74+0.50 


0.074+°:°°^ 


A1835 


Isoth. p 


O 1 q+0.09 

°--'-'-'-0.09 


0.543+°:°°} 


o 07+0.23 
O-<J'-0.23 


4 n4+o.08 


9 cc+O.ll 


0.158+°:°°^ 




Non-Isoth. p 


9.08+°!^ - 
61.6+1-^ ^ 


0.784+°:°}^ 




fi 40+0.30 
^•^^-0.33 


e- 77+0.70 


n 1 1 9+0.009 




100 kpc-cut 


33.6+i-° 


0.690+°:°°3^ 


11 4+°-'^ 


5 78+°-23 
'J- 'O-0.20 


5 Pifi+0-50 


n 104+0-005 
U.iU4_goo5 


^Inner 6 


c- 














'^Outer 9c. 















38 



Table 3. Sources of Uncertainty in /gas 
Measurements 



Source X-ray /gas SZE / 



gas 



Statistical Effects^ 

Kinetic SZE • • • ±4% 

Radio Point Sources • • • ±4% 

Asphericity ±20% ±10% 



Systematic Effects^ 



Instrument Calibration 


±6% 


±8% 


X-ray Background 


±2% 




HSE 


-10% 


-10% 


Isothermality^ 


-5% 


-10% 



^Statistical uncertainties average down by a fac- 
tor of v38 for the sample. 

^Systematic uncertainties do not average down. 
^Applies to isothermal 100 kpc-cut model only. 



Table 4. Model Parameters, Cluster Masses, and Gas Mass Fractions: Isothermal (3 model 

with 100 kpc cut 



Cluster 



(") 



ATo 
(mK) 



Tx 
(koV) 



xray 

AVs 
(IO^Mq) 



szc 

(IQI^^Mq) 



%tal 



xray 



J gas 



*cool/*Hubblo 



CL0016 

A68 

A267 

A370 

MS0451 

MC0647 

A586 

MC0744 

A611 

A665 

A697 

A773 

Zw3146 

MS1054 

MS1137 

MC1149 

A1413 

CL1226 

MC1311 

A1689 

RXJ1347 

MS1358 

A1835 

MC1423 

A1914 

A1995 

A2111 

A2163 

A2204 

A2218 

RXJ1716 

A2259 

A2261 

MS2053 

MC2129 

RXJ2129 

MC2214 

MC2228 



1.271 
1.60+ 



7+0.04 
-0.04 
,1+0.13 
-0.11 
2 04+"-^** 
73+002 
"■'•^-0.02 
, „c: + 0.11 
^■^'^-O.IO 

3.36+0-*3 



-0.36 
r,+0.12 

-0.10 
., O1+0.37 
■^■^-^-0.29 

"■'^-^-l.ll 

1+0.17 
-0.15 

'■°°-0.28 
3 + 0.02 

-0.02 
, + 0.11 

-0.08 
1+0.04 

-0.04 

4 31 + 0-59 
^■^-^-0.41 
R+0.76 

-0.48 
3 + 0.58 

-0.46 
R+0.29 

-0.26 
R+0.70 

-0.60 

l-'^-O.ll 

9 49+0-38 
^■^^-0.32 
o ,1+0.39 
■^■^-^-0.34 
,„ O + 0.40 
-^"■^-0.40 
, + 0.05 

-0.05 
3 + 0.10 

-0.08 
3 + 0.09 

-0.09 
7+0.35 

-0.22 
Q+0.03 

-0.02 
t-+0.22 

-0.13 
a + 0.19 

-0.14 
, + 0.96 

-0.77 
0+0.06 

-0.05 
a+0.25 

-0.20 
Q+1.10 

-0.83 
R+0.15 

-0.14 

2 36+0-1'' 
^■^"-0.13 



1.40^ 



2.31^ 



0.43: 
0.62^ 
0.81+ 



2.363 
1.03J 
6.261 
15.6+ 



1.623 
089J 
4.633 
1.371 
1.691 
0551 
1.483 
5.92I 
O.29I 
I.99I 
4.99I 
1.76+ 



42.81-.J 

rr 1+5.3 
OJ.l_5 

43.il=5;» 
Si-slal 

4o.oi^,;« 

24.1'' 



1+2.0 



,c 0+5.1 
''3-J_4.6 



25.8 
25.7 
17.7 
42.7 
38.7 
30.4 
57.7 
22.3 
44.7 
37.9 
16.2 
7.8 
48.4' 
17.7 
28.9 
33.6 
14.6 



+2.0 

2.0 
+ 2.1 
-2.2 
+ 2.9 

2.5 
+ 1.6 

1.5 
+ 2.5 

2.5 
+ 0.9 

0.9 
+ 8.8 

7.5 
+ 3.8 

3.4 
+ 2.6 

2.3 
+ 3.7 

4.3 
+ 3.8 

3.3 
+ 2.9 

1.5 
+ 2.1 

2.0 
+ 0.7 

0.8 
+ 1.8 

1.9 
+ 1.0 

1.0 
+ 1.2 

1.0 



44.oli;| 



57.5 
51.0 
57.6 
33.7 
70.3 
17.6 
59.1 
28.9 
27.0 
24.4 
25.4 



+2.7 

2.6 
+ 7.6 

6.6 
+ 1.3 

1.2 
+ 6.6 

5.7 
+ 1.7 

1.7 
+ 5.7 

4.5 
+ 9.1 
-8.4 
+ 3.5 

3.0 
+ 6.2 

4.8 
+ 2.9 

2.6 
+ 3.8 

3.3 



O.742IO 
0.76410 
0.720l° 
O.8I1I55 
0.773+° 



0.687' 
0.72310 
0.7231° 
0.632 



0.45410 
0.605155 

0.594155 
0.7361° 
0.8841° 
0.8771° 
0.6661° 
0.5361° 
0.7321° 
0.62ll° 
s+0 



0.688^ 



^-^-■^-2.8 
20.21^;^ 



0.65ll° 
0.6381° 
0.690+° 
0.6371° 
0.734+° 
1.0131° 
0.61ll° 
0.5381° 
0.6141° 
0.7651° 
0.6651° 
0.6631° 
0.6241° 
0.82ll° 
0.712!° 
0.61ll° 
0.710l° 
't°o 



0.522^ 



025 
.022 
046 
.042 
049 
.040 
039 
.037 
026 
.023 
023 
.020 
039 
.033 
028 
.027 
017 
.017 
006 
.005 

oil 

.011 

013 
.012 

009 
.009 

155 
.117 

097 
.078 

020 
.018 

015 
.015 

084 
.060 

030 
.022 

013 
.012 

007 
.007 

016 
.016 

007 
.008 

012 
.011 

009 
.010 

041 
.038 

044 
.036 

005 
.005 

047 
.034 

014 
.013 

097 
.064 

057 
.048 

023 
.018 

122 
.085 

040 
.033 

024 
.020 

038 
.034 

015 
.013 



1 0R1 +0.083 
1.0Dl_j, Qg3 

„ „i ,+0.099 
"•'-^^-0.103 

-0 692+"-°''^ 
U-D»^-0.080 
n SfiR+O-090 

-1 476+"-°^'5 
-^■^'"-0.086 
1 „i, +0.116 
l.olz_j, J23 

-0 651+0-°*! 
U-OO-1-0.084 

-1 299+0-134 
-0 761+0-°84 

-1 330+0-1" 
1 ,,„ + 0.125 

-1 101+"-"** 

-1 172+0-1*3 

l-l'^-0.144 

1 103+0.125 

-1-1O3-0.136 

-0 775+0-094 

"■' '^-0.102 

-1 166+0-112 

-1 037+0-150 



-1.7153 

-1.4891 
-1.651 + 



-0.150 
t;+0.202 
-0.234 
0+0.309 
-0.314 
1+0.130 
-0.137 
, „7„+0.294 
^■' ' ' -0.284 
„ 7,1+0.096 
-U- '^1-0.097 
1 ^0^+0-116 
-1.636_oii4 

-1 177+0.210 
^■^' ' -0.200 

-1 565+0-125 
1.0D0_j, J25 

-0 879+0-053 
"•^ '"-0.052 

-o.596lj;;ii^ 

-1 903+°-"^ 

-1 644+°-^'"' 

n s7n+0-079 

-O.S70_j, 075 

-0 656+0-1*8 
U-000-0.166 

-0 307+0.138 

1 17K+0-122 
-1-1'»-0.121 
-0 419+0-069 

-1 176+°-^2'' 
^■l'°-0.133 

-0 734+0-113 

-1 433+0-110 
l.«WJ_j, J26 
1 7,7+0.166 

"l-'^'-0.177 



10 46+°-01 
-^"■ID-0.60 

9 55+l-°0 
5.89l°„:f, 

s.erlS:!;, 

14 06+1-^8 
f. oc+0.46 
°--^^-0.39 
8 14+0-80 
°-l^-0.72 

6 79+0-41 
"■ '"-0.38 
„ 37+0.37 

"■•^'-o.ss 
10 21+0-''0 



8.163 
8.281 
9.77+ 



-0.65 

R + 0.56 

-0.52 

0+0.30 

-0.29 

7+1.10 

-0.94 

4 47+0.53 

*-*'-0.44 

9 93+0-80 
y.yo_Q gg 

7 51+0.35 
'■°^-0.29 
13 47+2-72 
lo.'l/_2 20 
_ jg + 1.53 
'■-'"-1. 16 

10 46+°-'"' 
-lU-''O-0.47 

sj+0.99 
-0.92 

3+0.88 
-0.74 

11 44+0-66 
i-1-11-0.56 

6 97+0-42 
°'"'-0.40 
9.48+0-35 



16.483 
8.93+ 



-0.29 
,+0.44 
-0.45 
(.+0.98 
-0.83 
1+0.39 
-0.38 
3+0.85 
-0.72 
n+0.41 
-0.37 
7+1.08 
-0.88 
1+0.44 
-0.36 
7 45+0-65 
'■*^-0.55 
4 77+0.71 
*■' '-0.58 
„ t;r+0.97 
°-^°-0.80 

6 69+°''^ 
°-°"-0.47 

10 22+0-99 

3+0.78 

-0.71 



8.223 
8.I5I 
I4.81I 
II.23I 
7.80I 
6.57I 
5.81 + 



a+0.29 

-0.29 
c^ + 0.34 

-0.31 
1+0.20 

-0.18 
s, + 0.21 

-0.20 
R + 0.31 

-0.30 

n+0.47 

-0.42 
2.26+0-13 



4.383 
3.651 
2.24I 
2.781 
4.761 
4.90+ 



3.07' 

2.36 

2.62 

4.39 

2.75 

4.42 

1.06 

1.19 

3.09 

2.63 

3. 

2.11 

5.06' 

8.84 

2.53' 

5.78' 

2.26' 



0.11 
+0.27 

0.25 
+ 0.11 

0.11 
+ 0.10 

0.09 
+ 0.28 

0.26 
+ 0.17 

0.15 
+ 0.11 

0.11 
+ 0.50 

0.49 
+ 0.14 

0.13 
+ 0.34 

0.30 
+ 0.11 

0.10 
+ 0.51 
-0.46 
+ 0.24 

0.21 
+ 0.16 

0.17 
+ 0.37 

0.35 
+ 0.19 

0.17 
+ 0.23 

0.20 
+ 0.10 

0.10 



4 85+0-12 
*-°°-o.io 



3.513 

2.19+ 



;.43"' 



1+0.14 
-0.14 
Q+0.27 
-0.23 
„ „o + 0.21 
»-"'*-0.20 
4 73+0-27 
*-''^-0.24 

3 00+0-13 
•^■""-0.12 
1 ,3+0.19 
-'■^'^-0.16 

1 S1+0-14 
1-81-0.14 

3 04+°-^'' 
•^■"^-0.19 

93+°-12 
O-"'5-0.12 

3 26+°-^° 
•^■^°-0.26 

2 58+°-10 
^■'^''-0.14 

3 95+°-33 
•3-»O-0.32 
, 01+0. 24 
^•81_o,22 



+ 0.26 
0.26 
+0.45 
0.44 
+ 0.32 
0.31 
+ 0.37 
-0.36 
+ 0.28 
0.28 
3 00 + 0.30 
^■^"-0.29 
2.49+0-32 



4.34 
3.07 
2.83 
3.51 
4.79 



-0.32 
K + 0.34 

-0.35 
3+0.26 

-0.26 
3+0.22 

-0.23 
+ 0.41 

-0.42 
1+0.39 

-0.39 
, + 0.47 
■^-0.47 

1 12 + °*^ 
l'l^-0.47 
3+0.26 
-0.25 
n + 0.30 
-0.30 
, Qo + 0.43 

, + 0.34 
-0.30 

3 + 0.38 
-0.38 

t; + 0.49 

-0.47 
, + 0.58 

-0.60 
3 + 0.29 

-0.29 
+ 0.36 

-0.36 
3 + 0.44 

-0.47 
5 11 + 0-12 
°'ll-0.42 

4 33+0-28 
4.JJ_g 27 
, 15 + 0.42 
^■l'^-0.42 

5 SQ + 0-55 
°'*y-0.53 
4 30 + °-^^ 

4 10 + 0-38 
^■10-0.38 
1 30 + 0.28 
-'-■•^^-0.28 

1 65+°-''l 
l''^^-0.58 
+ 0.40 
-0.38 
Q + 0.23 
-0.22 
+ 0.34 
-0.35 
7+0.33 
-0.33 
4 50 + °'^'^ 

1.0U_g 3g 

1+0.35 



3.363 
2.33I 
2.53I 
4.O9I 
3.641 
3.80I 



2.133 
3.00+ 



2.903 
2.33I 
5.851 
5.60I 
2.23I 
4.981 
2.63+ 



3.583 
I.39I 
3.281 
2.17+ 



3.843 



3.333 
4.32I 
2.O3I 

2.77I 
3.761 
5.98+ 



2.263 
2.14+ 



2.59^ 



3 + 0.39 

-0.37 
+ 0.91 

-0.74 
3 + 0.42 

-0.33 
7+0.37 

-0.33 
R + 0.53 

-0.46 
S+1.30 

-1.05 
2 48 + 0-34 
^■^"-0.28 
K + 0.42 

-0.35 
4+0.22 

-0.20 

2 00+0-14 
^■""-0.12 

3 47+0-40 
•^■^'-0.37 

Q+0.30 
-0.27 
3 62+0-22 
•^■°^-0.20 

74+0-45 
"■'^-0.38 

1 03+0-26 
l.uo_(j 20 
o 31+0.40 
^■'^-'--0.34 

2 15+01'' 
^■-'-°-0.15 
5 ,1+1.96 
'^■^-'--1.39 
O 17+0.74 
^■l'-0.53 
4gg+0.36 

n«+0-77 

8-08_o gg 

3 13+0.51 
•^■-'-^-0.41 

5 56+°™ 
''■°°-0.42 

1 94+0-20 
^■"^-0.18 

4 82+0-30 
*-°^-0.25 

^+0.50 
-0.48 

2 49+0-56 
^■^"-0.44 

5 49+0-24 
'^■^"-0.23 

4 96+°-80 
4-»o_0.64 

3 33+0.31 
ci.aa_g 27 
1 39+0-45 

1 79+0-27 
^- '"-0.24 

2 56+0-36 
^■^°-0.31 

1 ,3+0.38 
l-^'^-0.29 
o o,+0.57 

O n«+0-27 

3 85+0-68 
•^■'''^-0.59 

2 03+0-31 
^■"^-0.27 



4.74^ 



O.lSl"" 



131+0-007 

U.lal_g ggg 

n nR4+0-009 
0-084-0.008 
110+°-°" 

U.llU_g QU 

100+0-005 

U.lUU_g ggg 

n 1,7+0.009 
"■l^'-0.008 

n nR9+0-009 

091 + 0-007 

U.U»l_g gg7 

136+0-012 

U.iaD_g gjj 

110+0-006 

U.11U_Q QQg 

1+0.004 
-0.004 
127+0-007 
"■^^'-0.006 
106+0-006 
„ 1,,+0.004 
U-1^^-0.004 
144+0-014 
"■i*''-0.013 

115+0-014 
U.110_Q Q13 

134+0-008 

U.104_Q QQg 

„ i,,+0.004 
U--l^^-0.004 
075+0-015 
U-U'O-0.014 

097+°°^^ 
"■O"'-0.017 

102+°-°°* 
"■lO^_0.004 
109+0-006 

U.lU»_g QQ5 

n nsi +0-006 

0.0»l_g QQg 

104+"-°°^ 

U.lU4_g QQ5 

116+0-006 

U.llD_g QQg 

100+0-003 
U-iUU-0.004 
074+0-005 
"■O'*-0.005 

n nss+0-008 

0.0»»_g QQg 

147+0.003 
"■^^' -0.003 
096+°°°" 

U.U»D_g Q1Q 

0.0901°„:°q°q1 
0.0881°„:°q1? 
101+°°°'' 

U.lUl_g QQ7 

n iiq+0.008 

U.±±3_Q QQg 

076+0-012 

"■O'D_0.011 

116+°°" 
"■^^°-0.011 
124+0-008 
U-i^4-0.008 
103+0-009 

U.lUO_g QQg 

n 1 30+0.009 

U.100_g QQg 



130+0-016 

U.lCiU_Q QJg 

071+0-010 

U.U(l_g Qjg 

140+0-032 
"■^^"-0.029 
„ 1,7+0.020 
"■-'^'-0.018 
127+0-020 
"■-'^'-0.018 
055+0-014 

U.U»»_g QJ2 

100+0-010 

U.1UU_Q QJ7 

n 14S+0-032 
0-148_Q Q27 

109+°°" 
u.iuy_Q Q-^g 

126+0-014 

U.l^D_g QJ3 

„ j^_+0.019 
"■-'-'-'-0.017 
140+°-°^* 

U.l.iU_g Q21 

105+0-015 

U.1U0_Q QJ4 

153+0-034 

U.10Ci_Q Q2g 

207+0-054 
"■^O' -0.045 
130+0-022 

U.lJU_g Qjg 

139+0-022 
u.±oy_Q Q22 

056+0-024 

U.UOD_g Qjg 

107+0-047 

U.lU(_g Q35 

„ 11„+0.014 
0-118_o.Qi3 

069+°°" 

U.UD»_Q QJQ 

o.o7ii°:°Qf, 

0.0891°:°q1J 
„ 13t-+0.029 

U.lJO_g Q2g 

106+0-0" 

U.lUD_g Qjj 

092+0-013 

U.UMZ_g Qjj 

n no5+0.026 

0.0»5_g Q21 

107+0-011 
"■I'J'-O.Oll 

n OS7+0-014 
0.087_Q Q14 

123+0-0" 

U.lZO_g Qjg 

093+°°^* 

U.UMO_Q Q2g 

092+°°^'^ 

U.UMZ_g Q34 

+0.026 
0.023 

+0.039 
0.030 

+0.027 
0.023 



0.139 
0.113 
0.116 



104+0-022 

U.1U4_Q Qjg 

„ jj_+0.024 
"■^^' -0.020 

n 100+0.035 

0-l»»_o.031 



1.5 
1.4 
0.8 
2.1 
1.0 
1.0 
0.5 
0.3 
0.6 
1.1 
1.2 
1.0 
0.1 
2.8 
0.7 
1.9 
0.3 
0.9 
0.3 
0.3 
0.1 
0.2 
0.1 
0.1 
0.7 
1.2 
1.7 
1.2 
0.1 
1.2 
1.0 
0.8 
0.3 
0.9 
0.8 
0.1 
0.9 
0.9 



CO 
CO 



^Units are 10 erg cm arcmin s 



Table 5. Model Parameters, Cluster Masses, and Gas Mass Fractions: Non-isothermal 

Double-/? Model 



"■200 

(Mpc) 



(10 ^cm ••) 






ec2 
(") 



(1()1^ Mq) 



'"total 



xray 



CL0016 


2 12 + "-^'' 
^■1^-0.22 


A68 


2 {)2 + "'5'' 


A267 


1 r,r, + 0.25 
l'J'^-0.20 




1 si+f-i" 






MS0451 


1-^^-0.33 


MC0647 


2 02+"-^^ 


A586 


1 eo+^'i^ 

-^■""-0.15 


MC0744 


1 iq + 0-24 
1.19_Q j^3 


A611 


1 40 + 012 


A665 


2 11 + 0-24 


A697 


2 12+0-2S 


A773 


l.ob^O j^3 


Zw3146 


, „„ + 0.07 


MS1054 


1 94+037 


MS1137 


„ „„ + 0.15 
•J'^O-O.IO 


MC1149 


1 QK+Ol* 


A1413 


1 54+0.12 


CL1226 


, r,., + 0.34 


MC1311 


i.io Q j^3 


A1689 


1 «« + 0-07 


R1347 


2 -,,, + 0.07 
^■-^"■'-0.12 


MS1358 


, 00 + 0.15 

1.83 p j^g 


A1835 


1 Qe + 0.46 
1-^^^0.40 


MC1423 


1 30+'''"3 
-^■^"-0.04 


A1914 


1 89+"-^^ 


A1995 


2 B2 + 0-75 


A2111 


2 04+'''^3 
^■"*-0.44 


A2163 


«a+0.32 
2.89_o 23 


A2204 


2 43+0.21 
^■^■='-0.12 


A2218 


, 00 + 0.20 
1-83-0.15 


R1716 


, ,,+0.38 
-^■-^-^-0.16 


A2259 


1 47+0.23 
-^■^'-0.17 


A2261 


1 54+'''l'5 
-^■^^-0.14 


MS2053 


92+"'-'^ 


MC2129 


, 00 + 0. 17 
1-28_0.12 


R2129 


1 K8+0.14 


MC2214 


, „,+0.23 
1-O4_o.20 


MC2228 


1 94+029 
-^■^*-0.22 



7+0-57 

-0.43 
K + 2-09 

-1.67 
Q + 2-41 

-1.48 
,-f 1-98 

-1.40 
a + 1-66 

-1.07 
-,+1.78 

-1.35 
Q + 2.40 

-2.04 
r, + 3.81 

-2.47 
a + 1.18 

-1.01 
,r, + 0.76 

-0.64 
t; + 1.38 

-0.95 
« + l-41 

-1.18 
o (-0 + 0.20 
^■°^-0.19 
0.90+0-53 



1.27^^ 
2.961 
4.19J 
5.32+ 
2.181 
4.71I 
6.99I 
5.40I 
5.081 
2.35I 
2.75I 
3.78+ 



-0.30 

r,-f4.95 

-3.34 
K + 0.87 

-0.67 
+ 1. 87 

-1.80 
K+4.63 

-3.51 
,1+9.76 

-7.35 
c -^ + 1.07 
6.5610,73 
-,+0.59 

-0.42 
, + 0.69 

-0.53 
, + 0.13 

-0.09 

°'*'^-0.47 

ot. + 1.73 
8.85_^ 57 

1 70+f''5 
-^■'"-0.71 
, ,1+2.36 
^■^^-1.22 
-, t-7+0.56 
-^■^'-0.45 

e; + 0.33 
-0.44 

7+1.34 

-1.14 

K-r, + 6.76 

-4.36 



8.123 
2.461 
6.481 
5.761 
9.44+ 



6.413 
3.93I 
4.23+ 



4.25]; 
4.27I 
5.65+ 



3.74"t 



„ + 1.76 
-1.37 
6.86l?:f3 

6.5811,^06 

6 29 + '^'3^ 
°'^^-0.26 

4.17+ 



1.393 
O.82I 
1.04+ 



0+0.16 
-0.12 

,+0.03 
-0.03 

i+0.04 
-0.04 

54+"-''l 
"■^ -0.01 
-, „,.+0.03 
-^■•^^-0.03 

-^■'"-0.06 
-^ 7-^+0.05 



7+0.88 
-0.75 
e; + 2.62 
-1.69 
1 92 + + 05 
-^■"^-0.83 



4.95^ 



-0.05 

Q + 2.1 

-1.7 

a+1.00 

-0.61 

,+0.06 

-0.05 

a+0.02 

-0.02 

,+0.02 

-0.02 

7+0.5 

-0.5 

^+0.01 

-0.01 

n+0.09 

-0.09 

a+0.02 

-0.02 

0+0. 61 

-0.47 

2 49+0-1'' 

,+0.39 

-0.29 

,+0.11 

-0.10 

,+0.5 

-0.6 

,+0.7 

-0.5 

K+0.2 

-0.2 

7+0.4 

-0.3 

0+0.03 

-0.02 

n+0.01 

-0.01 

r+0.03 

-0.03 

i+0.01 

-0.01 

,+0.4 

-0.4 

,+0.01 

-0.01 

,+0.20 

-0.16 

^+0.04 

-0.04 

,+0.39 

-0.29 

^+0.14 

-0.11 

^+0.07 

-0.07 

,+1.6 

-1.2 

K+0.05 

-0.05 

-,+0-07 

'-0.06 



11.93 
3.281 
O.82I 
0.981 
O.93I 
15.73 
O.55I 
I.70I 
0.681 
4.781 
2.49I 
4.23I 
4.I3I 
25.33 
10. 2I 
12.61 
16. 7I 
I.49I 
O.90I 
O.55I 
O.94I 
21.13 
O.72I 
I.53I 
O.94I 
3.93I 
1.151 
I.75I 
14.13 
1.461 
I.31I 



,1+0.20 
-0.18 

7+0.18 
-0.15 

, + 0.18 
-0.16 

57+"-°'' 
"•'^'-o.oe 

, + 0.17 
-0.14 

n+0.14 

-0.13 

n+0.23 

-0.21 

n+2-7 

-2.2 

, + 0.94 

-0.66 

1+0.11 

-0.11 

92+0-1* 

,+0.18 

-0.15 

,+ 1.7 

-1.6 

, + 0.08 

-0.07 

^ + 0.32 

-0.32 

, + 0.06 

-0.06 

0+1.09 

-1.03 

+ 0. 42 

°-0.33 

n+1-20 

-0.95 

ir + 0.36 

-0.37 

0+1.9 

8.60l"l-i 

12.5101 

16.71=2-5 
-, ,7+0.11 
-'■''' -0.11 
1.09+ 



1.443 
O.77I 
1.22+ 



1.333 

i.iol 

I.50I 
12.03 
3.I2I 
O.91I 
O.92I 
I.31I 
14.23 
O.53I 
2.25I 
O.53I 
5.40I 

i.ssl 
4.701 
4.351 

16.9^ 



+ 3.2 
"-2.5 



Q+0.08 
-0.08 

57+0-1* 
"-="-0.12 
0+0. 09 
-0.11 



0.903 

,+ 1.9 
-1.7 
-,+0.10 
-0.09 
„;+0.42 
-0.37 
-. oQ+0.40 
l-^'^-0.34 
„+0.64 
-0.59 
77+0.60 
-"-0.45 

gg + 0.27 

-"^-0.23 

,+ 1.9 

-2.0 

sj+0.20 

-0.17 

1 58+0-28 
l.oo_o 22 



21.23 
O.81I 
I.35I 
I.39I 
4.14+ 



11.33 

1.58+ 



3.02 
4.59 
4.16 



+ 0.83 
0.67 

+ 2.89 
1.52 

+ 1.81 
23 



0.18 
0.19 



1.14 
1.20 



+ 1.5 



21.8 

, oo+o':i2 

•3-»8_o.l3 
3.30 
9.34 
4.08 

5.57 



+ 0.23 

0.24 
+ 0.21 

0.22 
+ 0.15 

0.14 
+ 6.60 

2.90 



-,+1.81 
1-1.15 
,+0.34 
^-0.30 



0.602"* 



4+0.44 
^-0.44 



759+0-034 
u.<o»_o 024 
707+0-031 
"■'"'-0.031 

694+0-029 
u.oy'i_o 029 
737+0.031 
"■''^'_o.027 
774+0.021 
"■' '^-0.019 

660+0-0" 
U.DDU_o 016 

627+0-015 
"-°^'-0.014 
584+0-072 
"-OO-*_0.049 
,+0.015 
-0.011 
661+0-01* 

U.DDl_o 014 

593+0-012 
u.aa3_o oil 

586+0-0115 

U.OOD_o 016 
n fiRR+0-006 
0.668_o 006 

1 743+0.188 
-L. ''±o_Q 236 

n fiKn+0-041 

0.680_o 039 
625+0-019 

U.D^O_o 017 

519+0-010 
u.oiy_Q 009 
704+0-044 
"-'"^-0.042 
615+0-025 
u.oio_Q 022 
0.889lO-0|3 

0.629lS-,Oog8 

676+0-010 
"-°'°-0.017 
„ 707+0.009 
"- '^' -0.013 
O.7O1IS-O- 

0.890lS:S? 
916+0-024 
u.yiu_Q 023 
596+0-032 

U.OMD_o 029 

559+0-006 
u.oo»_o 005 
741+0.028 
"-'*l-0.027 
746+0-015 

U.(4D_o 015 

540+0-052 
u.04U_o 043 

557+0-021 
"-^°'-0.019 
659+0-0*3 
U.DOM_o 032 

552+0-051 
u.ooz_o 055 

602+0-022 
"-o"^_0,020 
n 5S4+0-016 
U.584_o 015 

606+0-019 

U.DUDj, 019 
(1 "-,10 + 0.013 
O.ol8_j, 013 



n 4K-? + 0.046 
0.483_o 054 



930+0009 

U.»,3U_o 012 

570+0-083 
"-^'"-0.097 
077+0-073 
"■"' '-0.056 



n so, + 0.004 
O.S82_o 004 



752+0021 
"-'^^-0.021 



n 867+0-004 
O.S67_o 004 
942+0.003 
U.M4Z_o 004 
934+0-003 

940+0-001 

U.»4U_o 001 

976+0-002 
u.yio_Q 002 
009+0-021 
"-""»_0,008 



013+0-014 
U.U1,3_0 006 
957+0-003 
u.yo( -0.003 



0.034 
0.039 



0.009 
0.010 



48.1 
48.7 
40.6 
55.5 
34.0 
20.5 
31.9 
22.2 
23.3 
65.2 
42.2 
41.7 
26.0 
86.1 
14.6 
39.3 
36.6 
15.5 



+ 3.8 

2.6 
+ 3.3 

3.3 
+ 2.8 

2.7 
+ 2.8 

2.4 
+ 1.2 

1.2 
+ 1.0 

1.0 
+ 1.4 

1.4 
+ 7.1 

5.6 
+ 2.5 

1.6 
+ 2.0 

2.0 
+ 1.5 

1.5 
+ 1.9 

1.9 
+ 0.9 

0.9 
+ 6.1 

8.3 
+ 1.3 

1.2 
+ 2.1 

2.0 
+ 2.9 



(. + 1.6 
^-1.5 
9 2+7-2 

10418 

22.911;* 

37.411-8 

63.6ll-| 

36.01+;* 

67.7+ 



50.3 
49.9 
67.0 



11.2 
40.5 
44.6 
12.5 
18.6 
25.7 
22.3 
21.4 



+ 1.0 

1.0 
+ 1.5 

1.4 
+4.5 

4.3 
+ 1.2 

0.9 
+ 2.1 
-1.9 
+ 1.8 
-1.8 
+ 2.1 

1.7 
+ 2.8 

2.7 
+ 10.2 
-6.5 
+ 1.8 

2.1 
+ 1.2 

1.1 
+ 2.9 

2.6 
+ 1.3 

1.3 
+ 1.5 
-1.4 



+ 0.48 
0.56 
+ 0.50 
0.50 
+ 0.19 
0.19 
+ 0.19 
-0.19 
^ + 0.37 
-0.40 
4 90+0-39 
^-^"-0.36 

2 46+0-18 
^-*°-0.17 
+ 0. 24 
^-0.21 
^ + 0.12 
-0.11 
+ 0.12 
-0.12 
a + 0.35 
-0.38 
K + 0.18 
-0.18 
4 71+0.10 
^-'1-0.10 

1 06+0-25 

l.UD_o 25 
7 + 0.11 

-0.10 
, + 0.24 

-0.28 
t; + 0.10 

-0.18 
7+0.44 

-0.52 
, + 0.17 

-0.19 
K + 0.12 

-0.12 
e-+0.26 

-0.47 
, + 0.18 

-0.18 
+ 0.30 

-0.33 
K + 0.10 

-0.08 
, + 0.16 

-0.13 
+ 0.20 

-0.20 
1+0.35 

-0.37 
, + 0.49 

-0.54 
+ 0.29 

-0.22 
7+0.14 

-0.14 
, + 0.22 

-0.29 

t; + 0.13 

-0.13 
1+0.22 

-0.23 
+ 0.18 

-0.21 
o t., + 0.23 
a.oci_o 23 
+ 0.16 

-0.19 
^ + 0.37 

-0.44 

2 99+0-37 
^-^^-0.45 



4.63 
3.84 
2.35 
3.02 
5.25"* 



2.8 

2.553 

2.79I 

4.681 

3.06"* 



1.273 
3.O3I 
2.851 
4.I7I 
2.22I 
5.561 
9.55I 
2.82I 
6.49I 
2.561 
5.O2I 
3.79I 
2.31I 
8.62I 
5.O0I 
3.I7I 
I.33I 
I.95I 
3.34I 

o.gg"* 



2.803 

4.15"* 



K + 0.32 

-0.36 
, + 0.59 

-0.55 
,, + 0.35 

-0.34 
n+0.30 

-0.31 
n+0.43 

-0.49 
K + 0.41 

-0.38 
c: + 0.35 

-0.31 
7+0.30 

-0.31 
n+0.29 

-0.28 
+ 0.30 

-0.31 
a + 0.63 

-0.65 
„+0.50 

-0.46 
e; + 0.42 

-0.43 
, + 0.16 

-0.16 
Q+0.17 

-0.17 
K + 0.21 

-0.22 
+ 0.51 

-0.47 
c; + 0.33 

-0.32 
K+0.39 

-0.37 
^ + 0.46 

-0.49 
0+0. 60 

-0.58 
Q+0.30 

-0.28 
-,+0.43 

-0.41 
2 56+0-36 

4 29+0*2 
- ^^-0.38 
4 56+0-35 
''■OO-0.37 
2 39+0-*5 
^■■^^-0.45 
,1+0.66 
°-^l -0.72 
5.03 



4.763 
3.631 
2.761 
3.20I 
5.20I 
3.161 
2.I5I 
2.871 
2.40I 
3.O9I 
4.381 
4.34I 
4.25I 
I.03I 
1.691 
2.361 
3.181 
3.153 
2.461 
5.851 
6.40I 
2.39I 
6.41+ 



3.56 
1.16 
1.99 
3.50 
1.53 
3.41 
2.23 
4.49 
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Table 6. Model Parameters, Cluster Masses and Gas Mass Fractions: Isothermal (3 Model 

Constraints Derived from SZE Imaging Data 
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.^Q+0.71 
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^•' '-0.29 
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Table 7. Mean Gas Mass Fractions 



Method used in spatial analysis 



xX-ray 
J R 



'gas 



/■SZE 
J gas 



Joint fit to X-ray+SZE data, 100 kpc-cut isothermal f3 model 0.110 
Joint fit to X-ray+SZE data, non-isothermal double-/? model 0.119 
Spatial fit to SZE data only, isothermal f3 model 



+0.003 +0.006 
-0.003 -0.018 
+0.003 +0.007 
-0.003 -0.014 



0.116^ 



-0.005 +0.009 
-0.005 -0.026 
r. A 91 +0.005 +0.009 
U.iZi_QQQ5 _o.016 

n 1 9n+o-oo9 +0.009 

U.iZU_QQQg _o.027 
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Table 8. Comparison of Cool-Core and Non Cool-Core Results 



Analysis Method Dataset Cool-Core Non Cool-Core 

Joint fit X-ray+SZE 100 kpc-cut isothermal (3 model X-ray 0.110 ± 0.004 0.109 ± 0.004 

SZE 0.107 ±0.007 0.120 ±0.006 

Joint fit X-ray±SZE non-isothermal double-/? model X-ray 0.118 ± 0.005 0.120 ± 0.005 

SZE 0.113 ±0.007 0.125 ±0.007 

SZE-data only, 100 kpc-cut X-ray temperature SZE 0.098 ± 0.015 0.129 ± 0.011 

Joint fit X-ray±SZE no-cut isothermal p model X-ray 0.139 ± 0.009 0.106 ± 0.004 

SZE 0.158 ±0.014 0.118 ±0.006 

SZE-data only, no-cut X-ray temperature'^ SZE 0.148 ± 0.037 0.123 ± 0.009 

^ The large statistical uncertainties for this subsample reflect the small sample sizes: six clusters 
in the cool-core case and seventeen in the non cool-core case. 



Table 9. Constraints on rj, 



gas 



Method used in spatial analysis ^^s'^'^ Vgl^F 

Joint fit to X-ray±SZE data, 100 kpc-cut isothermal p model 0.65+°;°^ +°:°^ 0.65^°;°^ +o.o5 

Joint fit to X-ray±SZE data, non-isothermal double-/? model 0.71+^^ t^fg 0-^7 -om -om 

Spatial fit to SZE data only, isothermal (3 model ■ ■ ■ 0.64to:o3 -o:i4 
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Appendix A: Surface brightness profiles 
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Fig. 5. — : X-ray surface brightness radial profiles (points) with best fit 100 kpc-cut models 
(solid lines). Radii are calculated using a reasonable input cosmology of flM= 0.3, Qa= 0.7, 
and h = 0.7. The vertical dotted line is drawn at 100 kpc; data within this radius are not 
used in the fit. 
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Fie;. 5. — : Cont. 
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